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experimental work of others. Improvement in the condensation heat transfer coefficient was 
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tube cross-sectional area 
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condensate mass flux rate 
Nusselt number (aD,/4,) 
fluid pressure 
fluid pressure due to potential flow 
fluid pressure due to surface tension 
dimensionless parameter defined by Equation (II-85) 


radial dimension, cylindical coordinate system 
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wall 
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radius of curvature 

two-phase Reynolds number (p,U,D,/n)) 

vapor Reynolds number (p,U.D,/n,) 

tube wall temperature 

vapor Saturation temperature 

condensate film streamwise velocity 

vapor boundary-layer streamwise velocity 
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dimensionless vapor potential velocity defined by Equation (II-108) 
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condensate film normal velocity 
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dimensionless streamwise length defined by Equation (II-17) 
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dimensionless parameter defined by Equation (II-105) 
dimensionless parameter defined by Equation (II-34) 
dimensionless parameter defined by Equation (II-51) 

heat transfer coefficient 


vapor boundary-layer thickness 


A, vapor boundary-layer displacement thickness 

A, vapor boundary-layer momentum thickness 

§ condensate film thickness 

5° dimensionless condensate film thickness defined by Equation (II-64) 
n dynamic viscosity 

6 angular coordinate, cylindrical coordinate system 

K pressure gradient parameter defined by Equation (II-103) 
Kj suction parameter defined by Equation (II-103) 

ay, parameter defined by Equation (II-118) 

r thermal conductivity 

p fluid mass density 

o fluid surface tension 

v fluid shear stress 

Ts condensate/vapor interfacial shear stress 

Ts dimensionless interfacial shear stress defined by Equation (II-119,120) 
® parametric angle 

x(>) parametric radius defined by Equation (II-15) 
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I. INTRODUCTION 


A. MOTIVATION FOR INVESTIGATION OF ELLIPTICAL TUBES 

For many years, condenser design in the United States has been based on the Heat 
Exchange Institute (H.E.I.) or Tubular Exchange Manufacturers Association (T.E.M.A.) 
standards [1,2] which rely on tube-side empirical relations and averaged data to predict 
heat transfer parameters. Since there is a large degree of uncertainty in the accuracy 
of the prediction, the condensers are over designed to ensure sufficient margin of 
reliability. The result is an excessively large, expensive condenser for the desired 
thermal duty rating requirements. The majority of condensers currently used in the 
Navy were designed under these standards. 

Within the last two decades, computer modeling has been used to more accurately 
predict the heat-transfer coefficients of condensers. There are two basic approaches to 
computer modelling, an integrated and a differential approach, of which the former is 
most commonly used in condenser design (Walker [3]). Integrated approaches divide the 
condenser into zones with associated mean heat transfer and fluid properties which are 
typically determined from semi-empirical relations. These zones are integrated over the 
whole condenser and provide overall performance predictions with respect to shell side 
pressure drop and duty rating. Though easy to use, integrated methods are restrictive 
when used in conditions which are outside of the fluid parameter and geometry 
constraints of the empirical relations employed in the approach. Differential 
approaches solve the basic equations of fluid mechanics and heat transfer and provide 
an understanding of local variations within the condenser, but are too complex for use 


in condenser design. The improvements in accuracy as a rsult of computer design 


methods are a function of the empirical or theoretical models used in their development. 
Marto [4] provides a comprehensive overview of the state of computer design methods. 

The purpose of empirical and theoretical research into condenser performance is 
to more accurately predict condenser characteristics such that the designer can develop 
reliable, smaller, and hence less costly condensers which meet the same thermal duty 
requirements. Analyses of single tubes allow isolation of individual factors which 
affect condensation without the complexities which result from vapor flow over tube 
bundles or condensate inundation. Since the pioneering work of Nusselt [5], a 
significant amount of theoretical work has been done analyzing laminar film 
condensation on horizontal circular tubes (see Rose [6]), the results of which have 
greatly contributed to the ability to reasonably predict single tube heat transfer 
performance. 

Over the last two decades, major effort has been expended to study effects which 
enhance laminar film condensation heat transfer over that obtained from a plain 
circular tube. The majority of these techniques have focussed on controlling the 
thickness of the condensate film as this is the major resistance to heat transfer. The 
condensate film thickness generally increases with streamwise distance from the top of 
the tube and is dependent on both the rate of condensation and the interfacial shear 
between the vapor boundary-layer and condensate film. Additionally, at high vapor 
velocities, vapor boundary-layer separation occurs which results in a rapid thickening 
of the condensate film downstream of the separation point. Enhancement techniques 
include, but are not limited to, extended surfaces (fins), profiled tube surfaces (roped 
or corrugated tubes), and non-circular tube geometries. Marto [7], Bergles et al. [8] and 
Webb [9] provide comprehensive reviews on such enhancement techniques. These 
techniques enhance heat transfer through an increase in surface area to volume ratio 


and/or the use of surface tension to thin the condensate film. 


Non-circular geometries contribute toathinning of the condensate film by placing 
more of the surface in line with the direction of gravity; thus, making the average 
streamwise gravity component in the momentum balance of the condensate film larger 
than that obtained from a circular geometry. Dhir and Lienhard [10] applied a simple 
Nusselt type analysis to an arbitrary plane or axisymmetric body in which the 
streamwise gravity component varied with streamwise length. The results of this 
analysis gave expressions for condensate film thickness and Nusselt number which were 
identical to those derived by Nusselt except that the gravitational constant, g, is 


replaced by an effective constant, go, which is defined by 


x (g R)*3 
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where x is the streamwise distance from the leading edge of the surface and R is the 
radius of curvature of the axisymmetric body. These expressions can only be applied 
to systems for which curvilinear coordinates are applicable, i.e.,for bodies in which the 
radii of curvature are much greater than the film thickness. 

Shklover and co-workers [11,12] analyzed a horizontal cylinder with a logarithmic 
spiral type of surface curvature (see Figure I-1)such that more of the body surface is 
aligned in the direction of gravity and the radius of curvature is continuously 
increasing along the streamwise direction. They manufactured this tube by "mechanical 
deformation of a circular tube." The effect of this profile was to increase the 
streamwise effective gravity component over that found for a circular cylinders and 
to utilize surface tension effects on the pressure gradient to increase the mean film 
velocity, u(x,y),in the streamwise direction. Their analysis considered viscous, gravity 
and surface tension forces in the momentum balance of the condensate film and 


neglected vapor shear at the condensate/vapor interface. These effects resulted in a 





Figure [-1. Logarithmic Spiral Tube from Reference [12]. 


thinning of the condensate film and an increase of 20-30% in the overall heat-transfer 
coefficient when compared to a circular tube of the same surface area. These results 
agreed with experimental data. In analyzing the individual effects of gravity and 
surface tension, it was determined that surface tension was most significant over the 
top portion of the tube where dR/dx is largest and accounted for 10-15% of this 
increase. 

An elliptical tube with major axis aligned with the direction of gravity has a 
larger effective gravity than a circular cylinder, as well as an increasing radius of 
curvature over the top half of the tube. In addition, they are more practical from a 
manufacturing standpoint than the non-circular tubes analyzed above. Vapor 
boundary-layer separation should also occur at a point further downstream as compared 
to a circular tube and thus the rapid thickening of the condensate film would be 
delayed. Methods of manufacturing elliptical tubes, and associated problems will be 


discussed in Chapters IV and VI. 


Wallis [13] conducted flow visualizations (Figure [-2)of water around circular and 


elliptical tube bundles. These photographs provide a qualitative feel for the flow 
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Figure [-2. Flow Visualization, Water over Circular and Elliptical Tube Banks, 
from Reference [13}. 


characteristics around each type of tube bundle. In examining the photographs, it can 
be concluded that the separation point for the elliptical tube bank shifts downstream 
compared to the circular tube bank and results in smaller wake regions behind 
individual tubes. These effects should contribute to better heat transfer performance 
for the elliptical tube bank. Experimental studies by Joyner and Palmer [14] show that 
single phase flow resistance and pressure drop are significantly smaller for elliptical 
tube bundles as compared to circular tube bundles. In power condensers, shell-side 
(vapor) pressure drop has an effect on the overall efficiency of the power plant and on 


the thermal driving force for condensation (T,,,-T,,4)). Minimizing the vapor pressure 


cat 


drop improves power plant efficiency and maintains a relatively constant T,,, (for a 


large pressure drop, T,,, is reduced as vapor flows through the condenser resulting in 


Sat 
less thermal driving force). 

Studies involving single-phase heat exchangers and condenser-evaporators indicate 
that elliptical tube heat transfer performance is superior to a circular tube of 
comparable surface area. Ota and Nishiyama [15] conducted an experimental investi- 
gation of single phase forced convection heat transfer of air over an elliptical cylinder 
of minor-to-major axis ratio of 1:3 at various angles of attack. They concluded that 
elliptical cylinders gave improved heat transfer performance over circular cylinders at 
all angles of attack (as a result of fluid turbulence) in the range of Reynolds numbers 
studied (8000 - 79000). Moalem and Sideman [16] conducted a theoretical analysis of a 
horizontal elliptical tube in a condenser-evaporator used in desalinization plants. Their 
study showed a 10-20% enhancement for elliptical tubes as compared to circular tubes 
with the maximum enhancement achieved at a minor-to-major axis ratio (aspect ratio) 
of 1:4 where the major axis is aligned with the vertical, Huang and Mayinger [17] 
conducted an experimental free convection heat transfer study around elliptical tubes 
and found optimum improvement in heat exchanger performance for tubes with major 
axis aligned vertically and with an aspect ratio of 1:2. Merker and Bahr [18] used an 
analogy between momentum and heat transport to derive a semi-empirical relation from 
which to determine mass transfer rate from the fluid boundary-layer into the free 
stream. The results of their study showed improved heat transfer performance for heat 
exchangers in which the elliptical tubes were spaced wider in the longitudinal direction 
and more compact in the transverse direction. [In all these studies it appears that an 


elliptical geometry improves the heat transfer performance of a heat exchanger. 


Be 


SIGNIFICANT THEORETICAL STUDIES ON HORIZONTAL CIRCULAR 
TUBES 


The purpose of this survey is to provide a historical perspective to the 


understanding of laminar film condensation on a horizontal circular cylinder. Only 


those works which have made significant advances in this understanding are considered. 


The methods of some of these analyses will later be applied to formulate and solve the 


governing equations for the case of the horizontal elliptical cylinder. Figure I-3 


provides the geometrical layout and coordinate systems used in these analyses. 





Figure [-3. Geometry for Film Condensation on a 
Horizontal Circular Cylinder. 


Nusselt [5] used a simple momentum and energy balance to determine the heat 


transfer properties for condensation on flat plates and horizontal circular cylinders. 


[In simplifying these physical laws for solution, he made the following assumptions: 


(1) 
(2) 
(3) 
(4) 
(5) 


(6) 


The tube wall temperature is constant. 

The thermophysical properties of the fluid are constant. 

The temperature at the film/vapor interface is T,,.. 

The condensate film thickness is small compared to the radius of the tube. 
Heat transfer in the condensate film is one dimensional in the radial direction 
which implies a linear temperature distribution when the film is very thin. 
Convection in the condensate film is neglected. 

The only forces acting on the condensate element are due to viscosity and 
gravity. Inertial forces are neglected. 


(7) Flow in the condensate film is laminar. 
(8) The vapor is quiescent with no interfacial shear. 


The governing equations for this system become 


d*u 
"2 + p, g sind = 0 
(I-2) 
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where m is the condensate mass flux rate across the interface. The mean heat transfer 
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for a horizontal circular cylinder. 

Sparrow and Gregg [19,20] extended the Nusselt [5] analysis by applying boundary 
layer theory to the condensate film. This study incorporated convection and inertia in 
the energy and momentum balance while continuing to neglect interfacial shear. They 
determined that for practical engineering fluids, inertia and convection had negligible 
effect. 

Shekriladze and Gomelauri [21] included interfacial shear by considering 
momentum transfer across the interface due to the condensation process. They used an 
asymptotic expression for the interfacial shear based on an infinite condensation rate 
given by 

t, = m(U, - us) (Is) 
where Uj, is the streamwise velocity of the vapor at the outer edge of the boundary- 


layer. It was also assumed that Ug>»ug, such that the film velocity may be neglected and 


that inertia and convection effects in the condensate film could also be neglected. As 
noted by Rose [6], the validity of neglecting us requires that A,AT/n he, be small (large 
Prandtl number). For forced convection heat transfer on a horizontal circular tube, 


Shekriladze and Gomelauri determined the following expression: 


in = OS Re (I-5) 
where 
7 U_ D 
Re = Ara (I-6) 
nN; 


is the two-phase Reynolds number. When the effect of gravity cannot be neglected 
(mixed convection) the mean heat-transfer coefficient is given by 


Nu = 0.64Rey1 +/1 + 1.69 F (l=7) 


where 


nN, h, g D 
AMAT UZ 


F = (1-8) 


and measures the relative effects of vapor velocity and gravity. It should be noted that 
based on potential flow outside the vapor boundary layer, the interfacial shear will 
always be positive and therefore this method does not predict vapor separation as would 
occur in reality. Based on a minimum separation angle of 82° for single phase flow over 
a circular cylinder without suction and no heat transfer after separation, Shekriladze 
and Gomelauri conservatively estimated that the heat-transfer coefficient would be 
reduced by about 35%. Actual soma occurs at a point further downstream and 
results in a heat-transfer coefficient between the result for separation at 82° and no 
separation. These expressions have been found to be reasonably accurate for cases of 
high condensation rate and are useful for their simplicity. 


Fujii et al. [22] conducted an analysis of mixed convection condensation on 


horizontal cylinders using the two-phase boundary-layer equations. Their formulation 


neglected inertia, convection and pressure gradient effects in the condensate film. The 


governing equations for the model are given by 
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for the condensate film. The compatibility equations at the interface are given by 
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An approximate integral solution of the momentum equation for flow over a cylinder, 
with suction by Truckenbrodt [23], was modified to more closely agree with the 
numerical solutions of Terril [24]. This method enabled the determination of the 
interfacial shear and subsequent solution of the governing equations for the condensate 
film. This technique also predicted the point of vapor separation and subsequently, a 
more accurate heat-transfer coefficient. For locations downstream of the vapor 
boundary-layer separation point, the interfacial shear was assumed to be negligible and 
simple Nusselt theory used to determine the heat transfer. Local and mean Nusscit 
numbers were numerically determined for boundary conditions of uniform wall 


temperature and uniform wall heat flux. Their results were in fair agreement with 
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with experimental data. In addition to the dimensionless parameter, F, a new 


dimensionless parameter, G, was defined as 


q AT i, Py N; (I-10) 
nN, Ag \ Py 1, 


which is a measure of the condensation rate. As G-=, the results of Fujii et al. [22] 
agreed with the previous results of Shekriladze and Gomelauri [21]. Lee and Rose [25] 
noted an error in the formulation of Fujii et al. [22] in that the interfacial shear terms 
in the condensate film momentum equation should have been divided by 2. They re- 
analyzed the system of equations using both the modified and unmodified method of 
Truckenbrodt [23] and found little difference whether one used the modified or 
unmodified method of Truckenbrodt. 

To date, Gaddis [26] has conducted the most comprehensive study of condensation 
on horizontal circular tubes using the two-phase boundary-layer equations. His analysis 
neglected surface tension in the momentum equation as well as viscous dissipation and 
pressure in the energy equation of the condensate film. He determined that inertia 
effects were negligible for most media (with the exception of liquid metals) and that 
convection effects were somewhat significant for viscous liquids with Pr » 1, similar 
to the conclusions of the simpler analyses of Sparrow and Gregg [19, 20]. He identified 
three regimes of behavior for flowing vapors. For low vapor Reynolds numbers, Re,, 
gravity effects dominate the heat transfer and Nusselt’s analysis adequately represents 
the heat transfer behavior. For high Re,, condensation was shear controlled and the 
Nusselt number is proportional to the square root of Re, as found by Shekriladze and 
Gomelauri [21]. The intermediate region, in which gravity and vapor shear are both 
significant, results in no simple relation. Gaddis also analyzed several cases of flow 


separation (vapor boundary-layer and/or condensate film separation). For low Re, ( = 
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107), vapor boundary-layer separation occurred with no flow reversal in the liquid film. 
This result is expected as the pressure gradient effect on the condensate film, which 
tends to slow the condensate flow over the rear half of the tube, is a function of the 
vapor potential velocity. At moderate Re, (= 10°) vapor boundary-layer separation nd 
condensate flow reversal near the interface occurred with a resulting decrease in the 
mean fluid velocity and rapid thickening of the condensate film. At high Re, (= 10°), 
no vapor boundary-layer separation occurred but the condensate film separated at the 
wall causing a sharp reduction in mean condensate velocity with a rapid thickening of 
the film as a result of a significant, adverse pressure gradient. Fluid separation at their 
respective boundaries is defined as the condition at which the shear is less than or equal 
to zero. 

diMarzo and Casarella [27] conducted a similar analysis using a more general 
solution technique and arrived at the same results as Gaddis [26]. In their analysis of 
the flow separation phenomenon, they provided guidance on determining heat transfer 
performance once separation had occurred. At low Re,, gravity effects dominate and 
no flow reversals or thickening of the condensate film occurred. The use of a Nusselt 
type analysis with no interfacial shear would seem prudent in the region beyond vapor 
separation. At moderate to large Re,, pressure gradient effects dominate. The adverse 
pressure gradient over the back half of the tube causes flow reversal and a rapid 
thickening of the condensate film. In this case, it would be prudent to neglect heat 
transfer completely beyond the separation point of the vapor-boundary layer or 
condensate film. 

Rose [28] studied the effects of pressure gradient in the condensate film. In this 
case, the pressure gradient on the condensate film is due to the pressure gradient of the 


vapor, as determined by potential theory, which is impressed on the condensate film. 


MZ 


He simplified the analysis by using the Shekriladze-Gomelauri model [21] which ignores 
inertia and convection in the condensate film and uses the asymptotic expression for 
the interfacial shear. As a result of this analysis, a dimensionless parameter, P, was 
defined by 

_ Py he (I-11) 

p, A, AT 

which represents the pressure gradient effect (Gaddis [26] had a similar combination 
of dimensionless parameters which were equivalent to P). When P = 0, the governing 
equations reduce to those of Shekriladze and Gomelauri [21]. He concluded that the 
effect of including pressure gradient was to improve heat transfer over the forward 
half of the tube since the pressure gradient over this region is favorable (tends to 
increase the mean film velocity which reduces the condensate film thickness). As a 
result of the formulation, he found that for cases where P > F/8, a critical angle at 
some point on the rear of the tube was reached where d& /do ~~. In this case, it was 
not possible to obtain a solution over the entire tube. It was postulated that this critical 
angle might indicate some instability followed by some degree of waviness or 
turbulence. It was noted that this critical point was reached prior to the point at which 
the condensate film separated. For conditions which permitted solution over the entire 
tube (i.e. P < F/8), any increase in heat transfer achieved over the forward part of the 
tube was balanced by a decrease in heat transfer as a result of the adverse pressure 
gradient over the back half of the tube such that there was little change in the mean 
heat- transfer coefficient. Since the pressure gradient due to potential flow of vapor 


just outside the vapor boundary-layer is given by 


dU 
a Oe =a 
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which is symmetric around the surface of the tube from front to back, the favorable 
effect over the front half of the tube is exactly offset by the adverse effect over the 
back half of the tube. Referring to the studies conducted by Gaddis [26] and diMarzo- 
Casarella [27], in the region of moderate Re,, condensate film separation from the tube 
wall does not occur but the pressure gradient sufficiently retards the film such that 
flow reversals occur and the film still rapidly thickens. In the present study, the effects 
of pressure gradient using the analysis of Rose [28] resulted in a condensate velocity 
distribution which approached zero as the critical angle was reached. This phenomenon 
may be a result of the model used in formulating the problem. Since Shekriladze and 
Gomelauri [21] used an asymptotic expression for the interfacial shear, negative 
velocities cannot exist in the velocity profile. The positive velocity distribution is 
retarded by the adverse pressure gradient. In the case of P > F/8, the pressure gradient 
is large enough to cause complete stoppage of film flow at $, and the film thickness 
increases rapidly. 

Krupiczka {29] examined the effects of surface tension due to film curvature on 
condensation on circular cylinders. He used a simple Nusselt type model but included 
surface tension in the momentum equation. In his development he did not assume the 
film thickness was much less than the radius of the cylinder until after the inclusion 
of the surface tension term to account for the curvature of the film. The resulting 
equation was a second order ordinary differential equation which required two initial 
conditions. The first condition was given by the symmetry of the problem. However, 
the initial thickness was not obvious and was arbitrarily chosen to be that obtained 
from the Nusselt model. He concluded that the effect of surface tension was small on 
the forward part of the tube and increases in significance over the back part of the tube 


due to the rapidly changing film thickness. This significance was dependent on the 
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magnitude of dimensionless parameter, A, given by 


ana | 
Qa 
oa, (I-13) 


e; gr Pr 
which is dependent on the surface tension of the fluid through a "modified" Weber 
number (We = o/(pegr-)) and the radius of the cylinder, r. The conclusion was that 
surface tension can be significant for small diameter tubes and wires where, as r 
becomes smaller, the film thickness becomes relatively significant and hence the film 
curvature cannot be neglected. For most practical cases, r»6 and therefore the surface 
tension effect due to film curvature can be neglected. The effect of assuming the 
initial film thickness to be that given by Nusselt [5] was checked and verified to have 


minimal effect on the mean heat-transfer coefficient. 


C. SIGNIFICANT THEORETICAL STUDIES ON HORIZONTAL ELLIPTICAL 

TUBES 

Some theoretical work has been done on horizontal elliptical cylinders in a 
quiescent pure vapor using a Nusselt type model. The physical orientation and 
coordinate system is shown in Figure I-4. 

Cheng and Tao [30] approximated the surface of an ellipse by several circular arcs. 
They analyzed condensation on these arcs using the same assumptions as Nusselt. The 
ellipse was aligned such that the major axis, a, was aligned with the direction of 
gravity. From their numerical results they determined that the heat-transfer 
coefficient decreased with increasing eccentricity, k (defined as the ratio of the minor 
to major axis). Values of k were varied from that of a vertical flat plate (k = 0) to that 
for a horizontal circular cylinder (kK =1). The mean Nusselt number was determined by 


a surface area weighted average of the mean Nusselt number for each circular arc and 
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Figure I-4. Geometry for Film Condensation on a Horizontal Elliptical Cylinder. 
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C, is an integral constant which links each arc and / is the number of arcs used in the 
approximation. In the practical range of eccentricity (0.3 ~0.6),the mean heat transfer 
coefficient was increased by 10 to 18% over that of a circular tube with the same 
surface area. Ali and McDonald [31] conducted a similar type analysis as Cheng and 
Tao without the circular arc approximation as a first estimate for condensation on 
inclined circular tubes. 

Wang et al. (32}] used the Nusselt assumptions to analyze condensation on 
horizontal elliptical tube for which the major axis is oriented at an angle, a, witli 


respect to the vertical axis. They obtained an expression for the mean heat-transfer 
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coefficient given by 
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The results of this study showed that the maximum mean heat-transfer coefficient is 
obtained when the major axis is aligned with the vertical («@ = 0° ). These theoretical 
results were validated with experimental data for an elliptical tube with semi-minor 
axis of 8 mm and semi-major axis of 22 mm. Sheng and Cha’o [33] noted that the mean 
heat-transfer coefficient was incorrectly determined since the heat-transfer coefficient 
is averaged over the surface area and the radial distance, r, is not constant for an 
ellipse. In the course of this current study, the problem formulation of Wang et al. [32] 
was rerun with the correct expression for the mean heat-transfer coefficient (as later 
derived in Chapter II). The corrected mean heat-transfer coefficient of the elliptical 
tube analyzed by Wang et al. is 11.3% larger then a circular tube with the same surface 
area. Additionally, it should be noted that the heat transfer performance of this 
elliptical cylinder was better than the circular cylinder for angular orientations up to 
a = 50°. The corrected theoretical results are approximately 5% lower than those 
determined by Wang et al. and more closely agree with their experimental data. Figure 
(I-5) shows the heat transfer enhancement for elliptical tubes of varying eccentricity 
and orientation angle. 

Sheng and Cha’o [33] studied the effects of surface tension and variable wall 
temperature on condensation ona horizontal elliptical tube. The remaining assumptions 


of the analysis were the same as those above. For the wall temperature, a cosine 
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Figure (1-3). Elliptical Tube Mean Nu Compared to Circular Tube Mean Nu for 
varying Eccentricity and Orientation Angle. 


distribution of the form 

AT = (T,,-T,,(1 - A cos) (I- 18) 
was used (which has been shown (Memory and Rose [34]) to be in good agreement with 
experimental data) where the value of A depends on the ratio of the outside to inside 
heat-transfer coefficient. They determined that variable wall temperature affected 
local values but not the mean values of the heat-transfer coefficient. For the surface 
tension effect, it was assumed that the film thickness was much smaller than thc radius 
of curvature of the elliptical surface and thus the surface tension effect was due ely 
to the curvature of the tube wall. This is in contrast to the study of Krupiczka [| ‘or 
a circular tube where the surface tension was due to the curvature of the conde.. wte 


film. Over the back half of the tube the radius of curvature decreases with stream’. ise 
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distance resulting in a retarding effect on film flow. As the flow is slowed, the film 
thickness increases rapidly and at some critical angle becomes infinite in magnitude. 
This phenomenon is similar to that described by Rose [28] for the pressure gradient 
effect. The surface tension effect was determined to be significant for k < 0.6. In 
comparison toacircular tube, surface tension causes a favorable pressure gradient over 
the front half of the tube resulting in a greater film velocity, thinner film thickness and 
improved heat transfer. It has an opposite effect on the back half of the tube which 
tends to negates this improvement. The gravity component of the momentum balance 
is the driving force behind the enhancement in elliptical tubes. Consideration of 
surface tension results in a slight decrease in the mean Nusselt number as compared to 
the situation in which the surface tension is neglected. 

A potential advantage of an elliptical tube compared to a circular tube is the 
difference in vapor flow characteristics as a result of a better streamlined shape. 
Panday [35] developed an explicit numerical method for two dimensional film 
condensation and applied it to the case of downward flowing vapor over elliptical 
cylinders. Convection and inertia were included in the condensate film as well as 
surface tension and pressure gradient (as a result of potential flow of vapor outside the 
vapor boundary-layer). The interfacial shear was approximated using the asymptotic 
expression for infinite condensation rate. Several errors were found in the expressions 
for surface tension and pressure gradient as a result of an incorrect analysis of the 
differential streamwise length, dx. The first error is the result of assuming that dx = 
rdq@. This relationship assumes that the radial distance from the centroid of the ellipse 
is constant over the interval of the parametric angle. This fact is not true (as will be 
shown in Chapter II) and results in an error whose magnitude is dependent on the step 


size used in the numerical procedure. The second error involves the expression for r. 
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If the parametric angle is measured from the vertical axis of the ellipse (as appears to 


be the case in Panday’s expression for vapor velocity potential), then r should be given 


by 


x = b sind 


y = a cosh (I-18) 


r = yx? + y* = ajcos’h + ksin’h 


Panday used an expression for radial distance given by 

r= asin’ + k2cos’ (I-19) 
The magnitude of the resulting error is dependent on the eccentricity of the ellipse. 
Additionally, there appears to be an error in the gravity component of the momentum 
equation. Panday uses a body force given by pgsin®@, which is true only for circular 
geometries but is untrue for elliptical geometries. His results will therefore not account 
for the improved performance which are obtained by placing more of the surface in 
line with the direction of gravity. Panday’s conclusions were that the overall heat 
transfer was reduced for elliptical tubes at low velocities (as compared to a similar 
circular tube) due to a rapid thickening of the film at the stagnation point and the 
overall heat transfer is increased for high velocities due to increased interfacial shear. 
These conclusions are opposite to what one may expect. At low velocities, the film 
thickness near the stagnation point for elliptical cylinders is thinner than the equivalent 
circular tube as a result of increased effective gravity. This effect results in improved 
heat transfer. Pressure gradient effects due to potential flow and interfacial shear 
should be negligible since the velocity is small. At high velocities, the streamlined 
geometry of elliptical cylinders results in lower interfacial shear and pressure gradient 


effects as compared to circular cylinders. These effects should result in a thickening 
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of the condensate film as compared to circular cylinders and hence reduced heat 


transfer for the same effective diameter and free stream vapor velocity. 


D. OBJECTIVES OF THIS THEORETICAL STUDY 

The results of previous studies indicate that elliptical tubes can be used to 
increase condensation heat transfer in low vapor velocity condensers as compared to 
circular tubes with the same surface area. The present investigation examines the 
effects of vapor shear, pressure gradient and surface tension in laminar film 
condensation on a single horizontal elliptical tube with its major axis aligned with 
gravity and the free stream velocity. Interfacial shear is estimated using both the 
simple assumptions of Shekriladze and Gomelauri [21] and the more complex technique 
of Fujii et al. [22]. The latter case calculates the angle at which the vapor boundary- 
layer separates such that the effect of reduced drag on the mean _ heat-transfer 


coefficient can be evaluated. 
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Il. THEORETICAL DEVELOPMENT 

Following the theoretical developments used for the case of laminar film 
condensation ona horizontal circular tube, theoretical models are developed for laminar 
film condensation ona horizontal elliptical tube of eccentricity, k. This development 
starts with a Nusselt [5] type model in a quiescent vapor, and then adds forced and 
mixed convection as in the models of Shedkriladze and Gomelauri [21] and pressure 
gradient as considered by Rose [28]. The pressure gradient takes into account the 
effects of potential flow outside the vapor boundary-layer as well as surface tension. 
Finally, a model is developed which analyzes the vapor boundary-layer and boundary- 
layer separation following Fujii et al. [22]. Pressure gradient and surface tension are 
not considered in this model due to the complexities introduced by the use of the two- 
phase boundary-layer equation. Where possible these elliptical models are checked 
against existing theories for the "limiting" eccentricities of a circular cylinder (k=0), 


vertical flat plate (k=1) and horizontal flat plate (k~-). 


A. FACTORS RELATED TO ELLIPSE GEOMETRY 
Consider an elliptical cylinder whose cross-section is oriented such that the major- 
axis is aligned with the vertical as shown in Figure II-1. The eccentricity, k, of the 


ellipse is defined by 


(II- 1) 


Q le 


where a and b are the semi-major and semi-minor lengths, respectively. Functions 
related to the geometry of the ellipse are initially developed in a cartesian coordinate 
system (X,,y,) whose origin coincides with the centroid of the ellipse. This is then 


transformed into a cylindrical coordinate system (7,6) where r is the radial distance 
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from the ellipse centroid to a point on the ellipse surface and 9 is the angle measured 


from the vertical. The transformation equations are 


m7 sin8 


(II - 2) 
r cos6 


yy 


The elliptical surface can also be defined by a parametric angle, @, measured from 


the upper semi-major axis, such that the y, coordimate ona circle of radius a translates 





Figure II-1. Geometry for Film Condensation ona Horizontal Elliptical 
Cylinder with Major Axis Parallel to Gravity. 


to the y, coordinate on the ellipse and the x, coordinate ona circle of radius 5 translates 


to the x, coordinate on the ellipse. The transformation equations are given by 


x, = 5 sing 
(II-3) 


y, = 4 cosh 


Combining Equations (I-2) and (I-3) results tn a relationship between r,8@ and 4: 


sind = Fi sin8 (a) 


cosh = — cos6 (b) (il-4) 
a 
tang = a (c). 


It is assumed that the condensate thickness, 6,is much smaller than the radius of 
the elliptical surface. Therefore, the equations for the condensate film and vapor 
boundary-layer are developed using an orthogonal, curvilinear coordinate system, (x,y), 
where x is the streamwise distance along the elliptical surface and y is the distance 
normal to the surface. 


{; Radial Distance, r 


In cartesian coordinates, the surface of an ellipse is given by 


Z Z 
*, 71 08 (II-5) 
b2 a’ 


The radial distance, r, can be determined by 
2 2 4 
xy) 2¥xrtye Cr 


Using the parametric transformation equations (Equation (II-3)), an equivalent 


expression for the radial distance is given in terms of @ by: 


16) = af GS SG (1-7) 
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ie Radius of Curvature of an Ellipse, R 
The curvature, «(x,), and the radius of curvature, R(x,), of a function, f(x,), 


is given by 


Fe) | 
(1 + f@)yyP? | (II-8) 


R(x,) = 1/x(x,) 


K(x,) = 


where the prime denotes differentiation with respect to x,. Solving for y, as a function 
of x, using Equation (II-5), determining the radius of curvature using Equation (II-8) 
and transforming the results to parametric coordinates gives the radius of curvature as 
a function of @: 


R() = : [sin + k?cos*pP? . (II-9) 


a. Streamwise Length, x 
Since the radial distance to a point on the elliptical surface is not constant, 
the streamwise length is not proportional to 6 as it would be for a circle. Additionally, 
as eccentricity varies, the streamwise length for a given 6 is not the same. To obtain an 
expression for x in terms of parametric angle , consider first a point on the ellipse 
surface as shown in Figure II-2. Moving a small distance, dx, results in incremental 
changes, d@and dr. The resulting relationship between x,r and @ is given by: 
(dx? = (dr + (r dey. (II- 10) 
Expressions are now needed for dr and d@ in terms of d@. For dr, this is simply 
achieved by differentiating Equation (II-7): 


dp = 1_@ k=) sin2o i (I-11) 
2 Ioosth + Ksin2 


7095) 





Figure II-2. Enlarged View of Point on 
Ellipse. 


For d6@, taking Equation (II-4a), substituting Equation (II-7) for 7 and differentiating 


the resultant expression results in: 


cos8 d@ = ee , (I-12) 
[cos*> + k?sin*p)?” 


Combining Equations (II-4b) and (II-7), and substituting for cos@ in Equation (II-12) 
results in an expression for d@in terms of d¢@: 


ee (II- 13) 
cos?) + k*sin*> 


Substituting Equations (II-11) and (II-13) into Equation ({I-10) gives an expression for 


dx as a function of d¢: 


l one 
k*+—(k?-1)sin*2 ae 
4 ) vy) (IT- 14) 





do 


cos*) + k*sin*> 
This may be integrated to obtain the streamwise length. For compactness, a function 


7() is defined by 





pee “(k?-1)sin?26 (I-15) 


cos’) + k*sin’o 


In later developments, a characteristic length is used to non-dimensionalize 


the heat transfer parameters in the model. An effective diameter, D., is defined as the 


e? 


diameter of a circular cylinder having the same surface area as an elliptical cylinder 


and will be used as a characteristic length. The effective diameter is given by 


nD, = 2 [-x@ao (II- 16) 


Streamwise length is non-dimensionalized by: 


ox 
a = ; I-17 
= ( ) 


e 





This enables direct comparison between elliptical and circular cylinders since it 
represents an equivalent surface area. 
4. Streamwise Gravity Component, g, 
The component of gravity in the streamwise direction is tangent to the 
elliptical surface and is a function of streamwise location. A line tangent to the 


surface, defined by the slope of the surface, is obtained from Equation (II-5), 


dy, a x, 
i ne (II- 18) 
5 oP ne nh 
A unit vector in the streamwise direction is then given by 
: 1 [4 ies Qa 5 
Pe => _ 2 eee kay 
| "pe (II- 19) 





2  a* 2 
(b? - x;) + “a 


The streamwise gravity component is given by the dot product of the gravity vector and 


the tangent vector, 


ee 
\ 
OQ, 
b> 
ll 
1?) 
e/a 


Z (II-20) 
(b* - x) - ral 


Zt 


Transforming Equation (II-20) into parametric coordinates results in the expression: 


8, = & f,@) 

, (I-21) 

where fi(o) = ___ sing 
Vsin?d + k?cos7o 


Ss: Potential Flow at Ellipse Surface, U, 

When considering forced convection due to the flow of vapor over the 
condensate film, the velocity of the vapor influences the condensate film thickness 
through the vapor shear at the film/vapor interface. The vapor velocity at the interface 
will be dependent onthe vapor velocity outside the vapor boundary-layer. This velocity 
is determined from potential flow theory. Assuming that the film thickness, 6, and the 
vapor boundary-layer thickness, A, are much smaller than the radial distance, r, of the 
elliptical surface (and therefore may be neglected), the potential flow, Uy, about an 
elliptical surface with major axis aligned with the vapor free stream velocity, U,, is 
given by 


U, at a 


: V1 + kcot?o 


(11-22) 


Details of the derivation of this expression are provided in Appendix (A). 
6. Ellipse with Major Axis Perpendicular to Gravity 
For the case where the ellipse major axis is perpendicular to the vertical as 
shown in Figure II-3,the previously determined equation for streamwise component of 
gravity remains the same. For 0<k <1, the major axis is aligned with gravity. For 1 


<k<,the major axis is perpendicular to gravity. 
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Figure II-3. Geometry for Condensation on Horizontal Elliptical 
Cylinder with Major Axis Perpendicular to Gravity. 


B. FREE CONVECTION CONDENSATION ON AN ELLIPTICAL CYLINDER 
(NUSSELT [5] TYPE ANALYSIS) 


Consider a quiescent, saturated vapor at a temperature T condensing on a 


sat? 


horizontal elliptical cylinder of eccentricity, k,and semi-major axis, a, aligned with the 
direction of gravity as in Figure II-1. The same simplifying assumptions are made as 
proposed by Nusselt in his analysis of condensation on a horizontal circular cylinder. 
These assumptions are: 


(1) The vapor is pure and quiescent. 

(2) The tube wall temperature is uniform and constant. 

(3) The thermophysical properties of the fluid are constant and evaluated at some 
given reference temperature. 

(4) The temperature at the film/vapor interface is T,,,. 

(S) The condensate film thickness is small compared to the radius of the ellipse. 

(6) Heat transfer in the condensate film is one dimensional in the radial direction 
providing a linear temperature distribution across the film (based on 6«R). 

(7) The only forces acting on the condensate film are viscous and gravity forces. 

(8) Flow in the liquid film ts laminar with no waviness and no vapor boundary 
layer separation. 


A fluid element of unit width is analyzed in detail and shown in Figure II-4. 





Figure II-4. Condensate “Filme lement  fiormuleree 
Convection. 


Conservation of mass for the fluid element requires that, 


orf i dy 





de - pif. dy + 2 fo fu ‘) de+mde=0 | (I-23) 


where m is the local condensation mass flux rate. 


Defining the mean film velocity as 


Ll rs 
“= sl e dy, (11 -2em 


Equation (II-23) simplifies to 


u 


(#43) (II-25) 


Conservation of momentum for the film element is a balance of viscous and bodv 


forces which reduces to 


d’ 
n— + 8h) =O. (I-26) 
dy 


Cd 


Here, an assumption is made that the fluid is Newtonian, L.e., 


ay 3 (I-27) 


Conservation of energy is a balance of latent heat from condensation and heat 


conduction through the condensate film which reduces to 


A 
hm dx = Foss Ta dx 


(11-28) 
dy TygTra 
h,,6 
The boundary conditions for the momentum equation are 
Usa 
(I-29) 


Integrating the momentum equation subject to the boundary conditions results in 


a local film velocity of 


Fi) eS [sy : >| (1-30) 
ur 2 


and a mean film velocity of 


feels (I-31) 
1, ’ 


Combining the continuity equation, energy equation and expression for film mean 


velocity (Equations (II-25),(II-28) and (II-31)), results in 


4, AT alpr8 : 
ae - TEE (a) ?) (II-32) 
hy 5 dx| 3n, 


Using the relationship between dx and d@from Equations (II-14) and (II-15) and 


31 


carrying out the differentiation results in 





2 
1,?is iS 2 . (11-33) 
fe) + 3 A) 8 S 


4 
z - 3 (11-34) 
B, 
where B, is given by 
A,n, AT D 
1=— (II- 35) 
Pi hg 6 


Note that the equivalent parameter used by Nusselt [3] multiplies Equation (II-35) by 
a factor of three and uses radius in place of effective diameter. Substituting Equations 


(II-34) and (II-35) into (II-33) results in a first order ordinary differential equation for 








ine 
a, 4h) , _ 4x) (I-36) 
dp 3f(6) ' Df) 


This equation can be solved exactly by using an integrating factor which would 
require numerical integration or can be solved numerically using a forward stepping or 
propagation technique. The initial condition is determined based on the symmetry of 
the problem which requires that (dZ,/d)g.) = 0. Solving Equation (II-36) at ¢ = 0 


results in the initial condition 





3x(0) _ 3ak? 
Z, 90 7 TO = ge (11-37) 
D, f,(0) , 


a2 


The non-dimensional film thickness, 5, and film thickness, 6, are determined 
respectively by 
8°(b) = Zo)" and 
(II -38) 
8(o) = { Z,(b) B, } 


The local heat-transfer coefficient, a, and mean heat-transfer coefficient are given 


respectively by 


A, 
a(p) = —-~~ and 











5() 
ni II-39 
_ fra, 4 ae 
a = > ol = =e ald) x(>) db , 
0 AY 
and the local Nusselt number and mean Nusselt number by 
a(p)D D 
N. = e _ e 
u(p) 1, 3(6) 
(II- 40) 
a CDi rs 4 
Nu = pial hea 
u 1, x40 8(6) x(b) db 


For the case where k = 1 (circular tube); [,()=sin9g, f,’(@)=cos@, x(@)=a and D, = 
2a and Equation (II-36) reduces to 


a, Aco 7 _ 2 (11-41) 
; 1 : ? 
do 3 sind sind 





which is the same as that found by Nusselt [5] except for a difference in the definition 
of B,. 

The above development can also be applied to a horizontal elliptical tube whose 
Major axis is perpendicular to the direction of gravity as shown in Figure II-3. In this 


Case, D 1S ereater than a, i.e. k > 1. 
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CG FORCED CONVECTION CONDENSATION ON AN ELLIPTICAL CYLINDER 
(SHEKRILADZE-GOMELAURI [21] TYPE ANALYSIS) 


Consider a saturated vapor at temperature T flowing downward over a 


sat? 
horizontal elliptical cylinder with free stream vapor velocity, U,. The elliptical 
cylinder has eccentricity, k, and semi-major axis, a, parallel with the direction of 
gravity and vapor flow. The same assumptions as for the case of free convection are 
used here with the exceptions: 

(1) The force of gravity is neglected. 

(3) UZ > Uy 6. 


(4) The interfacial shear stress is approximated by an asymptotic expression 
assuming an infinite condensation rate, 


tT, = m(U,-u,) = mU, (II-42) 


Conservation of mass and energy are the same as previously derived (Equations 
(11-25) and (II-28)). The only forces acting on the fluid element are viscous forces. 


Thus, conservation of momentum reduces to 


2 
ny a a (II- 43) 
with the boundary conditions, 
Uy.g = 0 
(11-44) 


du 
nu (2) =t, =m U 
dy}, : 


Integrating the momentum equation subject to the boundary conditions results in a 


local velocity of 


m U 
y (I-45) 
Ni 





u(y) = 
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and a mean film velocity 





1D (II- 46) 
@ 


Using the expression for potential flow over an elliptical cylinder (Equation (II-22)) and 


defining a velocity function, f[,(), such that 


f(%) = colli _ : (II- 47) 


yl + kcot*o 


the mean film velocity can be expressed as 


U 
un = ae f(6) oa (11-48) 
ui 2 





Combining the continuity equation, energy equation and expression for mean film 


velocity (Equations (II-25),(II-28) and (II-48)) results in 


h, AT 1 _ d je 4, AT Us £06) é (II - 49) 
1, hy 2 





Using the relationship between dx and d@ from Equations (II-14) and (II-15), and 


carrying out the differentiation results in 


A, AT 1 p, aA, ATU. | dé 
en ee ‘ae 
fe 11 x() ) 
Defining a dimensionless parameter, Z,, as 
2 
(pe (I-51) 
B, 
where 
Began (I-52) 
Pp; U, 
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Equation (II-50) becomes 





dz, , fd) 4 xo) (I-53) 


—— Z, = 
db f,(o) D, f,(9) 





This problem is a first order ordinary differential equation which can be solved using 
a forward stepping numerical integration technique. The initial condition for the 
problem is based on symmetry which requires that (dZ,/d)j_9 = 0. Solving Equation 


(II-53) at @¢=0 results in the initial condition 


7 = 240) ime 
2 o=0 D, ££) D(1+k) 








(11-54) 


The non-dimensional film thickness, 5, and film thickness, 5, are determined 
respectively by 

8° = Z@) 

6 = /Z,() B, 


The heat-transfer coefficient for forced convection condensation can be considered in 


(II- 55) 


a different form. Combining Equations (II-40),(II-51), (I-52) and (II-55) results in 








Nu Re? = , and 
1) (I-56) 
Nu Re? = f° aad 
T™*9 8°() 


For the case where k = 1 (circular tube), f,(o) = 2sing, dx = ad and Equation (II- 


49) reduces to 


i eek Sn | | (II-57) 
a 


which is the same as the Shekriladze-Gomelauri [21] governing equation for a circular 


cylinder with no body forces. 
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D. MIXED CONVECTION CONDENSATION ON AN ELLIPTICAL CYLINDER 

(SHEKRILADZE-GOMELAURI [21] TYPE ANALYSIS) 

Consider a saturated vapor at temperature T,,,, flowing downward over a 
horizontal elliptical cylinder with free stream vapor velocity, U,. The elliptical 
cylinder has eccentricity, k,and semi-major axis, a, which is parallel with the direction 
of gravity and vapor flow. The forces acting ona film element are viscous and body 
forces. The interfacial shear is given by the Shekriladze-Gomelauri [21] model for 
infinite condensation rate as in the case of the forced convection model (Equation (II- 
42)). 

Conservation of mass and energy are as previously derived. Conservation of 


momentum for the condensate film reduces to 





d*u : 11-58 
Ni => + P18 Ad) = 0 (I-38) 
dy 
with boundary conditions, 
Uy. = 0 
G mu, (II-59) 
dy y=d Ni 


Integrating the momentum equation, subject to the boundary conditions, results 


in the local film velocity 








2 U 
u(y) = Ps 4) [5y-22) — = (0) y (II-60) 
N; 2 Ny 
and mean film velocity 
ec: Deel. Us 8 1-61 
u, - f,@) 3 = f,(%) a ( ) 


Where the expression for potential flow over an ellipse and its associated velocity 


function (Equations (II-22) and (II-47)) have been used. 
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Combining the continuity equation, energy equation and expression for mean film 


velocity (Equations (II-25), (II-28) and (II-61)) results in 


A, AT 83 A, ar U 
I —_ ie A) = oe, en = f,(0) ss (II-62) 





| pp aie neg 


Using the relationship between dx and dd@,and algebraically manipulating Equation (II- 


62) results in 








1 D 7 
¢ 6) 4 8 ({I-63) 
3° 2 x6) mak a ~ ae | 


where F is a dimensionless parameter relating free and forced convection and 6 isa 


dimensionless film thickness given respectively by 





7 De hy 8 
U22, AT 
(11-64) 
U : 
5° = 3 | re 
nN; D, D, 


Differentiating Equation (II-63) results in 


1 
So 2 ® 





|2F f0@ 6" oe 7 FAO) ~ + £9) & eo) 0| _ {ll-65) 


Regrouping the terms shows that the equation is a first order ordinary differential 


equation: 


{ 2F f,(o) 8° + £6) os are + |2F = - Ho) - 2X0) 1-66) 


é 


The initial condition for this problem is again based on symmetry, which requires that 
d8 /d=0 at 6=0. Inserting this condition into Equation (II-66) results in a fourth order 


expression for 6 





- F (0) 8" + f(0) 8” - 2 ” a“ (I-67) 


é€ 
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whose root of interest is given by 


1 
76) -f,(0) + (Go) + 4 2/3 F f,(0)} 7) a 


2 2/3 F f,(0)} 


1/2 
Bkl |(1+kY, 16 Fa _ +k 
4 k 3D, &k 


Equation (II-66) can be solved numerically to obtain the film thickness. The other heat 


[2 


(I-68) 





transfer parameters are determined as before (Equations (II-39) and (II-56)). 
For the case where F=0 and k=1 (forced convection over a circular tube), Equation 
(II-63) reduces to Equation (II-57) as expected. For the case where F - ~and k =1 (free 


convection over a circular tube), Equation (II-63) becomes 


1 D, d (2 ; 
2, = yee. SES a ae (11-69) 
5* 2x() _ we 


Upon substitution for the parameters F and 5° (Equation (II-64)), Equation (II-69) 


becomes 


2p ; 
_ 1 Pt hy 8 4) sing e : (I-70) 


1 
§ ada, AT n, a 


which is the same as Nusselt’s governing equation [5] for laminar film condensation on 


a horizontal circular cylinder. 


E. MIXED CONVECTION CONDENSATION ON AN ELLIPTICAL CYLINDER 
WITH SURFACE TENSION AND PRESSURE GRADIENT EFFECTS 
Consider a pure saturated vapor at temperature T,,,, flowing downward over a 
horizontal elliptical cylinder with free stream vapor velocity, U,. The elliptical 
cylinder has eccentricity, k, and semi-major axis, a, parallel with the direction of 


gravity and vapor flow. In addition to the assumptions used for mixed convection in 
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the last section, pressure effects in the condensate film are considered in the 
development of the governing equations. The pressure on a condensate film element 
will be a combination of that impressed upon the element by the potential flow of the 
vapor as well as that due to surface tension. A fluid element is shown in detail in 


Figure II-5. 





Figure [I-5. Condensate Film Element for Mixed 
Convection with Pressure Gradient. 


Conservation of mass and energy are as previously described. Conservation of 


momentum for the film element is a balance of viscous, pressure and body forces which 





reduces to 
2 
1, eet pO) — ee 0 (I-71) 
dy? dx 
with boundary conditions 
Lye = 0 
‘% m U, (II- 72) 
dy y=5 "1 
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As previously stated, the film pressure is a combination of the pressure due to 
potential flow, py, and the pressure due to surface tension, p,. Thus the total pressure 


gradient is given by 


dx dx dx 


dp _ Py , Po (II-73) 


The pressure gradient due to potential flow is given by 


wd = UP yan ay 2s (I-74) 
dx "9 de xh) db 


Substituting the expressions for potential flow and x() (Equation (II-22) and (II-15)), 


into Equation (II-74) results in 


U; 
ap, Py £0) (I-75) 


dx 2a 





where 





(1+k)* k? sin2o cos*h + k*sin*o 











f,(o) = (II- 76) 
[ sin’ + bcos’ ]* | p2 + 142-1) sin?2o 
4 
The pressure due to surface tension is given by 
oe ea?) 
" — -R() 
It is assumed that 5«R. Thus, the pressure gradient is given by 
ne ee a (11-78) 
ce: R? dx R? x(o) 4% 
Substituting Equations (II-9) and (II-15) into Equation (II-78) results in 
dp 30 
Somes Co (II-79) 
ile we 
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where 






cos*h + k*sin2p 


fae ; (k2-1)? sin?2p 


k (1-k*) sin2o 


[1-80 
[ sin’ + k2cos?o] *” cima 


Lb) = 


Integrating the momentum equation, subject to the boundary conditions results in 


a local film velocity of 


























2 
u(y) = _ 2 Hq, fi) + () + () (2, - z + -81) 
2 n, 
and a mean film velocity, 
: = Ato) + ( (| = + (6) 3. ({il-82) 
nN, @ 3 








Combining the continuity equation, energy equation and expression for mean film 


velocity (Equations (II-25), (11-28) and (II-82)), results in 














MAT 1d)! eis fio) * P1 Py Us 44 = SL" = 100) 
1 3 
he 8 dx 2a 1 3 oa 
_ pA ATU, 
Ti, ; 


Using the relationship between dx and dd,and algebraically manipulating Equation (II- 


83) gives 


| 














F f() + — * P Ab) + 3 (Plz f,(¢) 


1 D, 4! 
fay ~ X) db (11-84) 


1 A 
+ > f) 5°} 


where F and 8 are as previously defined in Equation (II-64). P and Bo are 
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dimensionless parameters defined by: 


ee altace 
p, A, AT 
(11-85) 
D2 
Roe P, 8 
0 


P is related to the pressure gradient and Bo, the Bond number, relates the inertia effects 
to surface tension effects. Differentiating Equation (II-84) results in the following first 
order ordinary differential equation: 


2 x) | 
D 


é 





D, D.) F 3 
er nae — Ese) tay a Bo 4) oe 


° D (I-86) 
+ £,(0) hile f F fio) + — P £@) 


+3 2. F Ko) ' + £@) | 

In the previous development, film separation and instability were not applicable 
since no effects were considered which could retard the film flow and (dé /dd)-@ only 
at ,=". However, in the present analysis, the decreasing potential flow and increasing 
radius of curvature of the elliptical cylinder causes the pressure gradient to have a 
retarding effect on the condensate film over the back half of the tube. At the point of 
condensate film separation from the tube wall, the film shear stress at the wall is 
negative which implies flow reversal (du/dy <0 at y =0). Differentiating Equation (II- 


81) and solving for the film separation criteria (du/dy = 0 at y = Q) results in the 


following condition for separation to occur: 


D p.y ; : 
2F f\(>) + = P f,() + 3 4 F 1) 6° +26) 8° <0 . (I-87) 
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Condensate film instability occurs when the rate of change of film thickness is infinite. 


From Equation (II-86), this condition occurs when 


D py : ; 
2F fd) + — P Ad) + 3 (4 & 1 a" + £(p) 8 = 0 2 


From a comparison of Equations (II-87) and (II-88), it is evident that the condition of 
film instability will occur first (has a smaller magnitude than the film separation 
criteria for all conditions) and is therefore the limiting condition. Rose [28] suggests 
that this instability may manifest itself in wavy or turbulent flow, or may signify the 
detachment of the liquid film from the tube surface. For the case where 1/Bo = 0 
(surface tension neglected), and using Equations (II-21),(II-47) and (II-76) for f,(), f(9) 
and f3() in Equation (II-88), it can be seen that the condition for which solutions can 


be determined over the entire tube surface occurs when 





P< ee F (II-89) 
For a circular tube (k=1), Equation (II-89) reduces to P< F/8, which was that found by 
Rose [28]. If the opposite of Equation (II-89) is true, then the pressure gradient effect 
is dominant enough to significantly retard the flow over the back half of the tube and 
result in a rapid thickening of the condensate film. Note that the condition for 
condensate film instability is a function of vapor velocity (F) since the pressure 
gradient due to potential flow is velocity dependent. 

For the case where P =0, and using Equation (II-88), the condition which allows 


solution of the momentum equation over the entire tube is now given by 





2 
AMigggie| GATE" ce (I-90) 
Bo 3\D,)} 1-k? 
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Note here that the condition for film instability is only a function of geometry. 
Therefore, for small values of eccentricity (approaching a vertical flat plate), the 
retarding effect of surface tension on the condensate film over the back half of the 
tube is significant and can result in a rapid thickening of the condensate film at @c¢. 

The initial condition for this problem is again based on symmetry, which requires 
that (dd /d)g-9 = 0. Combining this condition with Equation (II-86) results in a fourth 


order expression for 8: 


(ae, 2 ee a Can 2 2x0) | 1-91 
2 F FO) + —* P 0) (2) £ xo | 7 + 0) 8 * atte ) 


whose root of interest is determined by 
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2A 


6°(0) = 

Equation (II-86) can be solved numerically to obtain the condensate film thickness. The 
other heat transfer parameters are determined as before (Equations (II-39) and (II-56)). 
It can be shown that this formulation reduces to the models developed for circular 
tubes. Consider the case where k=1 and 1/Bo=0 (circular tube with no surface tension); 


Equation (II-84) reduces to: 


(11-93) 
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which is identical to that developed by Rose [28] where 


U,, P; 
2a n, 


- | ¢: (11-94) 


- 2 Rose 


6° = 6 








For the case where k=1, P=0 and 1/Bo=0 (circular tube, no surface tension and no 


pressure gradient), Equation (II-84) reduces to the Shekriladze-Gomelauri [21] model. 


F. MIXED CONVECTION CONDENSATION ON AN ELLIPTICAL CYLINDER 
WITH VAPOR BOUNDARY LAYER SEPARATION (FUJII ET AL. [22] TYPE 
ANALYSIS) 

The mixed convection model developed by Shekriladze and Gomelauri [21] has 
been used because of its simplicity and ease of solution. However, since the interfacial 
shear is approximated by an asymptotic expression based on an infinite condensation 
rate and uses potential flow outside the vapor boundary-layer (which is always 
positive), vapor boundary-layer separation does not occur. As described in Chapter I, 
Fujii et al. [22] modified the interfacial shear stress expression by simultaneously 
solving the boundary-layer equations for the condensate and vapor ensuring compatibil- 
ity at the condensate/ vapor interface. This technique allows for a more precise 
description of the condensation problem under forced convection conditions and 
includes vapor boundary-layer separation. 

Consider a pure saturated vapor at temperature T,,,, flowing downward over a 
horizontal elliptical cylinder with free stream velocity, U,. The elliptical cylinder is 
oriented such that its major axis is aligned with the direction of gravity and vapor flow. 
Due to the complexities of solving the two-phase boundary-layer equation, it is assumed 
that the only forces acting on the condensate film element are gravity and viscous 
forces. Surface tension and pressure gradient are neglected which is reasonable since 


they were found to have negligible influence under most conditions (discussed in 
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Chapter [V). Other assumptions are the same as in Section If.E. The vapor is assumed 
to be in thermal equilibrium such that the energy equation for the vapor need not be 
considered in the analysis. The only forces acting on an element in the vapor boundary- 
layer (see Figure II-6) are inertia, viscosity and pressure gradient, which is impressed 


upon the element by potential flow outside the vapor boundary-layer. The tangential 
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Figure II-6. Condensate Film and Vapor Elements for Mixed Convection using 
Fuji et al. [22] Type Model. 
velocity at the vapor/condensate interface is much smaller than the potential velocity 
outside the vapor boundary-layer and can therefore be assumed to be negligible. Based 


on the above assumptions, conservation of mass and momentum for the vapor boundary- 


layer are given by 
1S eget 
Ox dy 
dU 
pel ga . Cases Dow AU 
ax dy dx py, ay? 


0 (continuity) (a) 
(II-95) 


(momentum) (b) 
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For the condensate film, conservation of mass, momentum and energy are given by 


Ms aa + gfi(d) = 0 (momentum) (a) 
i (11-96) 
rm qe a Ay (Tso T att) 
m = at? u, \ = he (energy and mass) (b) 


where conservation of mass and energy have been combined as in previous sections. 


Boundary conditions for this system of equations are given by 


Uo = Vy» = 0 


(II-97) 
U ou 


= U ( | = 0 
=6+A o ? 
/ oy y=5+A 


Where A is the vapor boundary-layer thickness. The compatibility conditions at the 


interface require that 


Ou oU 
U =8 = u =8 = 0 ; Ny, ea = UE — a 


A, (T_.-T 
and pV g = sn vaw) : 
pie = h, 8 
fg 


(11-98) 


where the last compatibility relation is the condensation mass flux rate (i.e. for a given 
surface area, the amount of vapor condensed must equal the increase in the mass of the 
film). 

The technique for solving the vapor boundary-layer momentum equation involves 
an approximate integral technique developed by Truckenbrodt [23]. The initial problem 
development follows the method of Pohlhausen [36]. The continuity and momentum 
equations (Equations (II-95a) and (II-95b)), are integrated over the thickness of the 
boundary layer from y = 6 to y = +A and then combined to eliminate the normal 
component of vapor velocity, V, resulting in a momentum integral equation 

ld A, dU, 4 i 


— {u, a,} + + —* = + 


: = : (II-99) 
U; dx U, dx U, p, Us 
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where A, and A, are the vapor boundary layer displacement and momentum thicknesses 


given respectively by 


(II- 100) 


V,; is the normal component of vapor velocity at the vapor/condensate interface due to 
the condensation process itself and issometimes called the vapor boundary-layer suction 


velocity. Differentiating the A, term in Equation (II-99) and rearranging the equation 











results in: 
Gz, iene Us 2). H “Se (II- 101) 
dx U, dx : U, 0, ie 


With some further algebraic manipulation this yields 
2 
dx| n, U, 


To solve this expression, the vapor velocity distribution across the boundary-layer 


t, A, 


Ny U, 





2 
weer P| = (I-02) 
A, n, dx ny 








thickness must be determined. This distribution is a function of a pressure gradient 


parameter, x, and a suction parameter, x,, given by Tuckenbrodt [23] as 





Z 
n, a& Ny 


Truckenbrodt [23] defined a shear function f(x,x,) and a shape function H(x,x,) as: 





Gee A 
KxK,x,) = ee 7 and H(xx,) seq (II- 104) 


nN o A, 


and a dimensionless function, Z, given by 


(II-105) 
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Equation (II-102) becomes 


a 2U- 2 
Die a — ae ({I- 106) 


Based on existing theoretical data (Schlicting [37] and Torda [38]), Truckenbrodt [23] 


proposed a simple linear approximation for Equation (II-106) given by 


az _ YU, ’ 
‘a U, (0.441(1-2x,) - 6k]. (II- 107) 


Substituting Equations (II-103), (II-105) and (II-106) into Equation (II-107) and defining 


dimensionless velocities by 


a U - y 





U, = rr and V, = 7 (11-108) 
results in 
: D, dU 
ae) fo | 1 +2V,,/Re, | Se ae Z| : (11-109) 
dd gy. D x(b) dd 


The dimensionless suction velocity, V;s,is determined by transforming the compatibility 
relation involving the condensate mass flux rate in Equation (II-98) and is given by 


=y alee (II- 110) 


where Gis a dimensionless parameter which is proportional to the condensation rate and 


G = Bsa Frvau | Pr (11-111) 
n, hy py 7, 


Combining Equations (II-109) and (II-110), the governing equation for the vapor- 


is defined by 


boundary layer becomes 











x) | 0.441 | 1 - 2Ee | oe ee (11-112) 
; x(>) db 
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The initial condition for Equation (II-112) comes from symmetry and is given by 


ih) oe (II- 113) 
dd },.4 


Combining Equations (II-112) and (II-113) results in 





ss ES) (/Z@) )° + 0.882 G L@) YZ - 0441-0 | (II- 114) 
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This is a quadratic equation from which Z(Q) is determined to be, 
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=0 5*(0) e db 


The interfacial shear is determined from Equation (II-104). Truckenbrodt [23] 
provides a simple approximation for the shear function, f(x,x,), given by 


f (&x,) = 3.22 JK, (K,+K) (II- 116) 


where 


K, = 0.0682 + 0.174 x, . (II-117) 


x and x, are modified by combining Equations (II-103) and (II-105): 


K. = ?. wy Z and K, = Gag . (II- 118) 
x(b) db 5° 





Dimensionless interfacial shear is defined by: 


CESS (II-119) 


CC. 2. oo 
) 2 
p, UV. 


An expression for dimensionless shear can then be obtained by combining Equations (lII- 


104), (II-105) and (II-116): 


a U 
t, = 6 vk, (k,+k) (II-120) 
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The governing equations for the condensate film are solved in the usual manner. 
The momentum equation (Equation (II-96a)), is integrated to solve for condensate film 


velocity, 


: 28 4 sy : ar con (11-121) 
2 n 


The mean velocity is determined from Equation (II-24): 


: re aie oa (II- 122) 
27, 


Substituting Equation (II-122) into the energy/continuity equation (Equation (II-96b)) 


and multiplying both sides of the equation by 








pi Z 
1 Per (II-123) 
Pell ee) 
results in 
3 
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Carrying out the differentiation of Equation (II-124) and rearranging results in the 


ordinary differential equation, 





xo) _ F OF 1 Gi 50 
ds* _ dD, 3 4G & (II-125) 
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which, due to the symmetry of the problem, has the initial condition, 





Ca een (II- 126) 
do },4 
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Combining Equations (II-125) and (II-126) results in 








D 3 \db 4 G\ db 


7 


x0) _ F a) 60) - 7 E = 86°30) = 0 (II- 127) 
o=0 o=0 


from which 8 (0) can be determined by numerical methods. 

For the case of k = 1 (circular tube), Equation (II-125) reduces to the analysis of 
Fujii et al. [22] as corrected by Lee and Rose [25]. The differential equation for the 
vapor boundary-layer (Equation (II-112)) and for the condensate film (Equation (II- 
125)) are solved simultaneously using the compatibility relations of Equations (II-117), 
(11-118) and (II-120). At the point of vapor boundary-layer separation, the interfacial 
Shear stress becomes negative. Downstream of this point, it is assumed that the shear 
stress at the vapor/condensate interface is negligible and a simple Nusselt type analysis 
is used for the remainder of the elliptical surface. The heat transfer parameters are 


determined as before (Equations (II-39) and (II-56)). 
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lil, NUMERICAL METHOD 
In the theoretical development, the governing equations for each model were 
either reduced to a single first order ordinary differential equation (ODE) or to a 


simultaneous solution of a system of first order ODEs of the form 


D = faley) (III-1) 


Initial conditions were determined from the symmetry of the problem which required 
that both the slopes of the condensate film and vapor boundary-layer were zero at d= 
0. These systems of equations are solved using one of the forward stepping numerical 


methods. 


A. SOLUTION TECHNIQUE FOR THE NUSSELT [5] AND SHEKRILADZE- 

GOMELAURI [21] ANALYSIS METHODS 

Each analysis method resulted in a single first order ODE of the form given by 
Equation (III-1). The solution of this equation provides the condensate film thickness 
along the elliptical tube surface. Many numerical methods are available to solve this 
type of problem, each with its own advantages and disadvantages. For this study, a 
fifth order Adams method predictor-corrector algorithm was developed from Crandall 
[39]. The Adams methods are multi-step techniques which use finite difference type 
Operators incorporating previously determined points (hence the term multi-step) to 
determine the value corresponding to the next step. These methods provide the same 
accuracy and are more efficient than the Runge-Kutta (single step) methods but are 
restricted to a fixed step size. The Adams method employed in this study uses an 
Adams-Bashforth explicit method [39] to predict the value of the dependent variable at 


the next step and then an Adams- Moulton implicit method [39] to correct this predicted 
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value. This method is not self starting and requires determination of the first four 
points. To start-up the algorithm, lower order Adams methods were used to obtain these 


four points. The recurrence formulae for this algorithm are given by 


aw ae » B, fn(xpy,) (III -2) 


Where x is the independent variable, y is the dependent variable, h is the step size, s is 
the step number for which y is to be determined, k are the points determined in 


previous steps and 6, are the coefficients provided in Table III-1 and III-2. 
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TABLE III-1. RECURRENCE COEFFICIENTS FOR PREDICTOR FORMULA 


FORMULA 
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TABLE III-2. RECURRENCE COEFFICIENTS FOR CORRECTOR FORMULA 


=— a 
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(2) 
(3) 
(4) 













-1/12 8/12 5v2 O(h*) 


ae) | 1/24 -5/24 19/24 9/24 Oh?) 
~~" -19/720 106/720 -264/720 646/720 251/720 O(h*) 


The algorithm used in the solution for mixed convection condensation with surface 





tension and pressure gradient (Section II.E) is provided in Appendix (B). 
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For many applications, this type of start-up routine would normally not provide 
sufficient accuracy. However, no difficulties were encountered using this numerical 
procedure for solving the asymptotic interfacial shear stress models. Solution 
convergence was obtained for an angle step size of 0.1°. The difference between an 


angle step size of 1.0° and 0.1° was less than 1%. 


B. SOLUTION TECHNIQUE FOR THE FUJII ET AL. [22] ANALYSIS METHOD 

Reduction of the governing equations resulted in a system of two simultaneous 
ODEs, one for the momentum thickness of the vapor boundary-layer and one for the 
condensate film thickness. The Adams method previously described did not provide 
sufficient accuracy to solve this two-phase boundary layer problem. In particular, the 
system of equations exhibited a stiffness problem, 1.e.,the numerical method could not 
accurately solve the vapor boundary-layer ODE over the first several steps even with 
an angle step size of 0.001°. To obtain sufficient accuracy, a much smaller step size 
would be required which would significantly reduce the efficiency of the computer 
algorithm. Lee [40] indicated that similar difficulties had been encountered when Lee 
and Rose [25] solved the same equations for a horizontal circular cylinder. Lee and 
Rose [25] used a Runge-Kutta method which permitted variable step size. Thus, the step 
size could be reduced sufficiently to obtain the required accuracy for the vapor 
boundary-layer ODE over the front of the tube and then increased over the remainder 
of the tube to improve the efficiency of the algorithm. Utilizing this numerical, 
procedure, they were able to solve the system of equations over most of the range of 
dimensionless parameters (F and G). However, in some cases (specifically large values 
of F and G), the system of equations were still too stiff for solution. 

To overcome the problem of stiffness, a numerical problem solver (IVPAG) from 


the International Mathematical and Statistical Library (IMSL) [41] was used. This 
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algorithm uses Gear’s stiff method [42] to efficiently solve systems in which stiffness 
is a problem. However, as in the study of Lee and Rose [25], this method could still not 
solve the two-phase boundary-layer equations over the entire range of dimensionless 
parameters considered, but did cover the practical ranges of the parameters. 

It was originally conceived that the computer algorithm would progress around 
the tube surface until the criteria for vapor boundary-layer separation was reached 
(dts/d@ < 0). At this point, the model would shift from solution of the two-phase 
boundary-layer equations to the simple Nusselt model. However, the numerical solver 
encountered difficulties in solving the system of equations as it approached the point 
of vapor boundary-layer separation. Therefore, the solution method was changed toa 
two step procedure. In the first step, the system of equations were solved until the 
computer algorithm indicated a problem had been reached (usually the algorithm 
indicated a stiffness problem had been encountered). The interfacial shear stress data 
was analyzed to verify that this problem was the result of rapidly changing conditions 
associated with vapor boundary-layer separation. Therefore, to facilitate the solution, 
vapor boundary-layer separation was assumed to have occurred if the interfacial shear 
stress at the last point obtainable was less than ten percent of the maximum interfacial 
Shear stress observed over the surface. At this point, the problem solver was within one 
to two step sizes of the actual separation point (for a step size of 0.1°, the computed 
separation point was therefore within 0.2° of the actual separation point). The second 
step involved re-running the analysis using the computed separation point calculated 
in the first step as the point at which the two-phase model was switched to the simple 
Nusselt model. A converged solution was obtained for an angle step size of 0.1°. The 


algorithm for solving the two-phase model is provided in Appendix (C). 
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IV. RESULTS AND DISCUSSION 
When comparing solutions for an elliptical tube with those of a circular tube, 
equivalent surface areas have been used. The effective diameter of the equivalent 
circular cylinder is given by Equation (II-16). The dimensionless streamwise length x , 
defined as the ratio of the streamwise length, x, to the half perimeter length and given 
by 


—— (IV -1) 





is used when comparing local values of film thickness and heat-transfer coefficients. 
Its use enables direct comparison between elliptical and circular tubes since it represents 
an equivalent fraction of the total perimeter from the top of the tube to the 
dimensionless point. 

A practical range of eccentricities for an elliptical tube with major axis aligned 
vertically is 0.3<k <0.6. The lower limit of 0.3 is based on discussions with a tube 
manufacturer (Reference [43]) in which the elliptical tube is formed by pressing a 
circular tube with a roller assembly to obtain the proper major and minor axis 
dimensions. Roller contact points would be evenly spaced to maintain the elliptical 
curvature and to prevent an overly flat surface. The upper limit is chosen as the 
eccentricity in which noticeable effects were observed in the heat transfer 


characteristics (as discussed in this chapter). 


A. EFFECT OF GRAVITY 
As previously noted (Section I.A),some improvement in condensation heat transfer 
is expected for an elliptical tube in quiescent vapor compared to a circular one since 


more of the surface is aligned with gravity. By increasing the "effect" of gravity on the 
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condensate film, the mean condensate velocity at a given streamwise location is 
increased which results in a thinning of the condensate film and increased heat 


transfer. The gravity function, [,(), defined in Equation (II-21) is shown in Figure IV- 
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Figure I1V-1. Gravity Function, f[,() vs. Streamwise Length, x . 


As the elliptical tube approaches a flat plate configuration (k -0), the effect of gravity 
increases to sweep the condensate film along the tube surface. The effect of gravity on 
the local film thickness is shown in Figure IV-2 for two eccentricities. As can be seen, 
the condensate film is thinned over the front and rear portion of the tube compared to 
a circular tube but is slightly thicker in the middle region. The thickness of the film 
is controlled by the rate of condensation (which is dependent on condensate film 
thickness) and the condensate velocity (which is also dependent on condensate film 


thickness and gravity). Near the top of the elliptical tube, the larger gravity component 


increases the condensate velocity relative to that of a circular tube resulting in a 
thinner condensate film. This thinner film, however, results in an _ increased 
condensation rate which tends to thicken the film further downstream. Since the 


elliptical tube has a greater condensation rate near the top of the tube, the film thickens 
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Figure [V-2. Local Film Thickness for Varying Eccentricity (k). 


more rapidly compared to that ona circular tube causing the relatively thicker film in 
the middle region as seen in the figure. At 90°, the gravity component is the same as on 
a circular tube and so has no relative thinning effect. Over the rear half of the tube, 
the reduced condensation rate (as aresult of this thicker condensate film) and increased 
velocity due to the larger gravity component relative to a circular tube results 1 
continued thickening of the film on the elliptical tube, but at a slower rate than that 
on a circular tube. Eventually, the elliptical tube film thickness is again thinner than 


on a circular tube for the same dimensionless streamwise length. The condensate film 
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thickness is therefore determined by a balance of the effects of two effects: increased 
condensate velocity (which tends to thin the film) and increased condensation rate 
(which tends to thicken the film). 

The overall effect on the mean heat-transfer coefficient (given as a mean Nusselt 
number) is shown in Figure IV-3. To compare the elliptical tube with a circular tube, 


the leading coefficient of the Nusselt solution (Equation I-3b) is plotted against 





Figure [V-3. Mean Nu Leading Coefficient for Varying k. 


eccentricity. Solutions for a circular cylinder (k = 1.0) and a vertical flat plate (k - 0, 
L = 2a) agree well with the results of Nusselt [5]. For the horizontal flat plate (k - »), 
the heat-transfer coefficient approaches zero as would be predicted by the Nusselt 
model. For the practical range of eccentricities (0.3 < k < 0.6), it can be seen that the 
effect of placing more of the tube surface in the direction of gravity is to increase the 


mean heat-transfer coefficient by 7% for k=0.6 and 13% for k = 0.3. Conversely, by 
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placing more of the tube surface perpendicular to the direction of gravity (i.e., 
approaching a horizontal flat plate), the mean heat-transfer coefficient decreases as 


expected. These results are in agreement with earlier elliptical tube studies [32,33]. 


B. EFFECT OF VAPOR VELOCITY 
ie Asymptotic Interfacial Shear Stress Approximation 
The streamlined shape of the elliptical tube has an effect on the vapor flow 
over the tube. The vapor velocity function, f,(), defined by Equation (II-47) is shown 
in Figure [V-4 versus dimensionless streamwise length x. As seen in the figure, the 
elliptical tube experiences higher vapor velocities than a circular tube at the front and 


rear of the tube but a lower vapor velocity in the central region. Using a balance of the 
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Figure [V-4. Vapor Velocity Function, f,(), for Varying k. 
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factors which thicken and thin the condensate film. for the case of pure forced 
convection (F -0,no gravity effect), the higher relative vapor velocity at the top of the 
elliptical tube results in a larger interfacial shear stress and hence a thinner condensate 
film and larger condensation rate compared to a circular tube. With more condensate 
flowing into the middle section of the elliptical tube and a lower relative vapor 
velocity, the condensate film thickens more rapidly than in the case of acircular tube. 
As the condensate flows over the rear of the tube, the decreased condensation rate (due 
to the thicker film in the middle region) and increased vapor shear results in thinner 


film compared to the circular tube. Figure IV-5 shows these effects in detail for a 
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Figure IV-5. Effect of Vapor Velocity on Local Film Thickness for a Circular 
Tube and Elliptical Tube (k = 0.6). 


circular tube (kK = 1) and an elliptical tube (k = 0.6) using the asymptotic interfacial 


shear stress (Shekriladze-Gomelauri) approximation and the two-phase boundary layer 
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Shear stress (Fujii et al.) approximation (to be discussed later). Whercas the gravity 
effect (f,()) described in Section IV.A was always larger for the elliptical tube (see 
Figure IV-1), the vapor shear effect f,() is larger or smaller depending on the section 
of the tube to be analvzed (see Figure 1V-4). The overall effect of vapor shear using ‘he 
Shekriladze-Gomelauri [21] model is to reduce (by about 2% in the practical range of 
eccentricities) the elliptical tube mean heat-transfer coefficient relative to an 
equivalent circular tube. 

In the mixed convection region, the mean heat-transfer coefficient of an elliptical 
tube is increased or decreased compared toa circular tube depending on the relative 
magnitudes of the vapor shear and gravity effects. This relative magnitude 1s measured 
by dimensionless parameter, F. For large F, gravity is dominant and the heat-transfer 
coefficient trend is described in Section IV.A . For small F, vapor shear is dominant 


and the trends described in this section determine the heat transfer. Figure IV-6 shows 
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Figure IV-6. Effects of Vapor Velocity on Mean Nu Using the Asymptotic 
Interfacial Shear Stress Approximation Analysis Method. 
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the effect of F on the mean heat-transfer coefficient for varying eccentricity. As F - 
eo, corresponding to free convection condensation, solutions for elliptical and circular 
tubes match those reported in Section IV.A. For the case of a circular tube for all 
values of F,the present numerical solutions are within 0.4% of an equation by Rose [28]: 


(1+344F24 Fy 


Nu Re™'? = 
Diz Two-Phase Boundary-Layer Shear Stress Approximation 

Vapor boundary-layer separation is not predicted by the asymptotic shear 

stress approximation since the interfacial shear is based on potential flow outside the 
vapor boundary-layer which is always positive. Solution of the two-phase boundary- 
layer equations for the condensate and vapor, with matched shear stress at the interface, 
allows determination of the vapor boundary-layer separation point and its effect on the 
condensation heat transfer. The analysis also introduces a dimensionless parameter G 
which is proportional to the condensation rate. Figure IV-5 compares the asymptotic 
shear stress approximation with the two-phase boundary layer shear stress approxima- 
tion for a large condensation rate (large G) for a circular and elliptical tube (kK = 0.6). 
Over the forward and middle sections of the tube, the general trends are as described 
in Section IV.B.1 for the asymptotic shear stress approximation. Over the rear of the 
tube, however, as the separation point is approached, the condensate film thickens more 
rapidly due to the reduced vapor shear effect (t - 0). It can be seen that for the 
elliptical tube, the separation point occurs further downstream than for the circular 
tube and thus this rapid thickening is delayed. By delaying separation, more of the tube 
surface has a thinner film. This serves to increase the heat-transfer coefficient for an 


elliptical tube compared to a circular tube if it were the only effect considered. 
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In total, however, the mean heat-transfer coefficient is affected by a 
combination of gravity, vapor shear and boundary-layer separation. Figure IV-7shows 
the effect of F, G and k on the vapor boundary-layer separation point (ay As 
previously explained (Chapter III), solutions could not be obtained for all possible F and 
G combinations (particularly large F and G). In general, as F and G increase (less vapor 
shear and higher condensation’ rate), or k decreases (more elliptical), the vapor 
separation point shifts toward the rear of the tube. For low condensation rates (small 
G) corresponding to low vapor suction, vapor separation occurs near the positions 
obtained for single phase separation without suction. Little change occurs in the 
separation point location as vapor velocity is varied. For large condensation rates (large 
G), the vapor boundary-layer experiences high suction which shifts the separation point 
downstream (i.e., aS G increases, xo increases). As vapor velocity is decreased 
(increasing F), the separation point shifts further downstream. As k decreases (more 
elliptical), the vapor boundary-layer separation point moves further downstream when 
compared toa circular tube under all conditions of F and G. At high F and G, the vapor 
boundary-layer separation points come together at the rear of the tube for both the 
elliptical and circular tube. 

Figure 1V-8shows the influence that the above effects (vapor shear, gravity, 
condensation rate and eccentricity) have on the mean heat-transfer coefficient. The 
solutions for k =1 agree with the results of Fujii et al. [22] as corrected by Lee and Rose 
[25]. Similar trends were obtained for the elliptical tubes. For large values of 
condensation rate (large G), the two-phase boundary-layer shear stress analysis agrees 
closely with the asymptotic shear stress analysis method as expected since the latter 
assumes an infinite condensation rate. For decreasing k, the results for an elliptical 


tube tend to show a reduction in heat transfer at low F (compared to a circular tube) 
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VAPOR SEPARATION POINT 
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Figure [V-7. Vapor Separation Point Location with Varying F, G and k. 
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Figure IV-8. Effects of Vapor Velocity and Vapor Separation on Mean Nu using 
the Two-Phase Boundary-Layer Approximation Analysis Method. 


and an increase in heat transfer at high F, as previously mentioned. The cross-over 
point from reduction to improvement occurs at increasing values of F as G decreases. 
For large values of F, the results converge to the Nusselt analysis. This shift in the 
cross-over point is attributed to the change in the vapor boundary-layer separation point 
with varying condensation and vapor velocity (varying Gand F). As seen in Figure IV- 
7, for large condensation rate (G=5.5) and low vapor velocity (large F), the large vapor 
suction causes a delay in vapor boundary-layer separation such that it occurs at 
essentially the same streamwise location for all eccentricities. As velocity increases (F 
becomes smaller) the separation point shifts forward at different rates dependent on 
eccentricity which is evidenced by the diverging curves in Figure IV-7c. Since early 
vapor boundary-layer separation results in a reduction in heat transfer due to rapid 


thickening of the condensate film, the circular tube (k = 1) experiences more of a 
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reduction in heat transfer compared to an elliptical tube. The effect of delaying 
separation is to cause a thinning in the mean film thickness for the tube. In the 
balance, where the streamwise component of gravity and the shift in vapor boundary- 
layer separation are positive effects (i.e., both contribute to a relative thinning of the 
elliptical tube condensate film compared to a circular tube) and vapor shear is a 
negative effect (i.e., lower shear contributes to a relative thickening of the elliptical 
tube condensate film compared toa circular tube), delays in vapor separation move the 
cross-over point from gravity dominant to shear dominant heat transfer. For the case 
of low Gin Figure I[V-7a, there is no noticeable shifting of the vapor boundary-layer 
separation point as F is varied. Thus the vapor boundary-layer separation point 
contribution does not change with vapor velocity (F) to the balance and the cross-over 
point occurs at lower velocities (larger F). The overall effect of vapor shear using this 
analysis is typically a slight reduction (< 2%) in the heat-transfer coefficient for an 
elliptical tube compared to an equivalent surface area circular tube, dependent on the 
magnitude of the condensation rate parameter, G. 

Figures IV-9 and IV-10 show the effects of vapor boundary-layer separation 
on the local film thickness and heat-transfer coefficient for a given dimensionless 
streamwise distance. For low F and G, the separation point occurs relatively early, but 
is significantly delayed when using an elliptical tube. As F increases (vapor velocity 
decreases), the separation points remain at their same respective streamwise locations. 
The areas under the curves in Figure IV-10 represent the total heat transfer. Bearing 
this in mind, the comparison between elliptical and circular tubes is made a little easier. 
At high F (for given G), the area under the curve for an elliptical tube exceeds that for 
a circular tube while at low F, the opposite is true. However, for a given eccentricity, 


the values of oe are the same for varying F (as noted earlier). 
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As G increases (increasing condensation rate and suction), vapor boundary-layer 
separation is further delayed (compared to low G) and is virtually eliminated at high 


F. Separation is always delayed for an elliptical tube compared to a circular tube. 
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LOCAL FILM THICKNESS 
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Figure [V-9a and b. Dimensionless Local Film Thickness for Varying G, F and k 
from the Two-Phase Boundary Layer Shear Stress Method of Analysis. 
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Figure IV-9c and d. Dimensionless Local Film Thickness for Varying G, F and k 
from the Two-Phase Boundary Layer Shear Stress Method of Analysis. 
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Figure 10a and b. Local Nu for Varying G, F and k from the Two-Phase Boundary 
Layer Shear Stress Analysis Method. 
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Figure 10 cand d. Local Nu for Varying G, F and k from the Two-Phase Boundary 


Layer Shear Stress Analysis Method. 
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C. EFFECT OF SURFACE TENSION 

As noted in Section II.E, a condition exists for which d& /d@ becomes infinite at 
some critical angle @-. What is actually happening can be discerned from Figure IV-11 
which shows the surface tension function f,() versus x for varying eccentricity. 
Surface tension causes a favorable pressure gradient over the front half of the elliptical 
tube and an adverse pressure gradient over the back half. The severity of the pressure 


gradient is localized to small regions at the top and bottom of the elliptical tube where 
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Figure I[V-11. Surface Tension Function, f,(@) for Varying Eccentricity. 


change in surface curvature is most severe. Therefore, it is most significant for small 
values of eccentricity. The above analysis explains the reason for condensate flow 


instability discussed in Section II.E and predicted by Equation (II-90). 


(s 


Further discussion of surface tension is restricted to those cases in which a 
solution could be obtained over at least 90% of the tube surface. This restriction was 
arbitrarily chosen based on the small amount of heat transfer which occurs on the lower 
part of the tube and would therefore have minimal impact on the accuracy of the 
results. Figure IV-12 shows the effects of surface tension on the mean heat-transfer 
coefficient for k = 0.6 and 0.4. The Bond number, Bo, gives the relative effect of 
inertia to surface tension. Values of 1/Bo of 0.01 (typical for steam) and 0.001 (typical 
for refrigerants) as well as 1/Bo = 0 have been shown on the figure. A number of 
interesting points can be highlighted. Firstly, the surface tension effect is much smaller 
for highly wetting refrigerants, as expected. Secondly, inclusion of surface tension 
leads to a small decrease in the mean heat-transfer coefficient (< 2%) over the whole 
range of F for the practical range of eccentricities, suggesting that any thinning of the 
condensate film over the top half of the tube is more than offset by a thickening over 
the lower half. As eccentricity decreases, this discrepancy is accentuated. Finally, 
surface tension effects are felt more in the free convection region (high F) than in the 
forced convection region (low F) as film thickness becomes dominated by vapor shear. 

Krupiczka [29] used the curvature of the condensate film surface to analyze 
surface tension on a horizontal circular cylinder. He determined that the effect of 
surface tension was only significant over the bottom portion of the tube. This fact 
implies that since the film is very thin over the top of the tube, its curvature is very 
close to that of the tube surface itself, which, for a circular tube results in no surface 
tension effect (constant radius of curvature). As the condensate flows around the lower 
half of the circular tube, the film curvature no longer follows the tube surface. Rather, 
the radius of curvature is increasing resulting in a small improvement in the heat- 


transfer coefficient. A similar analysis can be applied to the elliptical tube. Over the 
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Figure [V-12. Effect of Surtace Tension on the Mean Nu for & = 0.4 and 0.6. 
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top half of the tube where the condensate film is thinnest, modeling the surface tension 
effect using the curvature of the tube surface is valid. However, the adverse pressure 
gradient over the lower half of the tube is probably not as severe as seen in Figure IV- 
11. Over this portion of the tube, the film is thicker and its decreasing surface 
curvature is not as severe as the tube surface curvature. The conclusion drawn from 
this analysis is that surface tension does not have a significant effect on the mean heat- 


transfer coefficient for an elliptical tube. 


D. EFFECT OF PRESSURE GRADIENT 

The effect of pressure gradient was analyzed using the asymptotic shear stress 
approximation (similar to Rose [28] for a circular tube). In this analysis, the pressure 
gradient in the condensate film was assumed to be due to the pressure gradient 
impressed on the condensate by vapor potential flow. Asin the case of surface tension, 
there isarelationship between parameters P and F (Equation II-89) for which a solution 
could not be obtained for the entire tube. Figure IV-13 shows the pressure gradient 
function versus x for varying k which sets up a favorable pressure gradient over the 
top of the tube and an adverse pressure gradient over the lower half of the tube. The 
effect of eccentricity is to shift the point of maximum pressure gradient to the front 
and rear of the tube. Thus, @c shifts downstream for an elliptical tube compared to a 
circular tube for similar P and F. This is a result of the more streamlined shape and the 
smaller resultant pressure drop over the streamwise length of an elliptical tube. 

As in the analysis of surface tension, discussion is limited to those cases for which 
solutions could be obtained over the entire tube. When pressure gradient was included 
in the momentum equation of the condensate film, solutions agreed with the results of 
Rose [28] for a circular tube, yielding a maximum 5% decrease in the heat-transfer 


coefficient over the entire range of F compared to the case with pressure gradient 
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Figure IV-13. Pressure Gradient Function, f3($), for Varying Eccentricity. 


neglected. Figure IV-14 shows the mean heat-transfer coefficient for a circular tube 
(k = 1) for varying F and P. The favorable pressure gradient increases the mean 
condensate velocity resulting in a thinner film and an increased mean heat-transfer 
coefficient over the top of the tube. This increase is offset by the adverse pressure 
gradient which cad the condensate film velocity over the lower half of the tube, 
thickens the film and decreases the mean heat-transfer coefficient. For an eccentricity 
of 0.6, Figure IV-15 compares the effect of pressure gradient for an elliptical tube with 
those for acircular tube for varying F and P. For both the circular and elliptical tubes, 
the effect of pressure gradient provides aslight reduction in condensation heat transfer. 
Figure [V-1S5indicates that the reduction in heat transfer was less for the elliptical tube 
compared to the circular tube as a result of the more favorable potential velocity 


distribution around the ellipse. 
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Figure [V-14. Effect of Pressure Gradient on the Mean Nu for a Circular Tube 
(k=1). 
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Figure [V-15. Effect of Pressure Gradient on Mean Nu, & = 0.6, Relative to a 
Circular Cylinder. 
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E. EFFECT OF VAPOR PRESSURE DROP 

The previous results indicate that vapor shear decreases the heat-transfer 
coefficient for an elliptical tube when compared to a circular tube. However, a 
consideration not yet discussed is that the vapor pressure drop across an elliptical tube 
is significantly smaller than a circular tube for the same free-stream velocity. 
Alternatively, for a given pressure drop, the free-stream velocity is greater for the 
elliptical tube. Since very little data is available for two-phase drag coefficients around 
tubes, single-phase data was used where the drag coefficient for the tube is due entirely 
to form drag. Values of Cp for a circular (kK = 1) and an elliptical tube (k = 0.5) were 
estimated to be approximately 1.2 and 0.6 respectively (White [44]). For steam with T,,,= 


60°C condensing on a horizontal tube with T 40°C and U,= 25 m/s for the circular 


wall 
tube, the corresponding vapor velocity for the elliptical tube, for the same pressure 
drop, was calculated to be U,= 35 m/s. The resulting values of F are therefore 0.0257 
and 0.0132 for the circular and elliptical tubes respectively. Using the asymptotic shear 
stress assumption for both values of F gives an increase in the heat transfer for an 
elliptical tube of 16.3% when compared to a circular tube. Using the two-phase 
boundary-layer shear stress assumption gives a corresponding increase of 17.1%. 
Therefore, though the effect of vapor shear alone results in a reduction of the heat- 
transfer coefficient (as discussed in Section IV.B), when the streamlined shape of the 
elliptical tube and its effect on the vapor pressure drop is taken into account (allowing 


for higher U, for a given pressure drop), there is an increase in the heat-transfer 


coefficient when compared to a circular cylinder. 
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F. INSIDE HEAT TRANSFER COEFFICIENT 

So far only the outside heat-transfer coefficient has been discussed. In real 
condensers, it is the overall heat-transfer coefficient which is important and so 
consideration must also be given to the single-phase convective heat transfer occurring 
inside the elliptical tube. A preliminary survey of the literature (Incropera and Dewitt 
[45]) shows that for turbulent flow conditions where the Prandtl number is greater than 
0.5 (which is true for most practical condensers), correlations developed for circular 
tubes may be used with good accuracy to approximate the inside heat-transfer 
coefficient for an elliptical tube if the hydraulic diameter is used in place of the 
circular diameter. The hydraulic diameter is given by 


_ 4c (IV-3) 


5 





D,, 


where Ac is the cross-sectional area of the tube and 1s the inside perimeter of the tube. 
The perimeter for an elliptical tube and an equivalent surface area circular tube is the 
same. The cross-sectional area of an elliptical tube is given by 

ef (IV-4) 
and is smaller than an equivalent surface area circular tube. Inside heat transfer 


correlations for turbulent flow heat transfer inside a circular tube are typically of the 


8 
Nu | = wt x Re®®| = acepliash sia 
A n 


Thus, the inside heat-transfer coefficient is related to the hydraulic diameter by 


a & De : (IV -6) 


form [45] 


(IV-5) 
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As an example, consider an elliptical tube with k = 0.5 and a major axis with a = 0.01 
m. The equivalent circular tube would have a diameter of 0.0154 m. The ratio of the 


heat-transfer coefficient of the elliptical tube to the circular tube is given by 














-0.2 -0.2 
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@ circle (D 2)’ 








For this simple analysis, it can be seen that the inside heat-transfer coefficient of the 
elliptical tube is 3.5% greater than the equivalent surface area circular tube. This 
result, when added to the increase in the outside heat-transfer coefficient (for a given 


pressure drop), would increase the overall heat-transfer coefficient of the condenser. 
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V. CONCLUSIONS 


Analyses of laminar film condensation on a horizontal elliptical tube were 
conducted under conditions of free and forced convection, the latter using both an 
asymptotic interfacial shear stress approximation as well as solution of the two-phase 
boundary-layer equations. For the asymptotic shear stress approximation, the effects 
of surface tension and pressure gradient were also included. Where possible, these 
analyses have been validated against existing solutions for laminar film condensation 
on a horizontal circular tube and a vertical flat plate. 

Whether or not condensation heat transfer is improved for an elliptical tube 
compared to an equivalent surface area circular tube is dependent on a balance of the 
effects of gravity, vapor shear and vapor boundary-layer separation and the influence 
these have on the thickness of the condensate film. Under quiescent vapor conditions 
(no vapor shear or boundary-layer separation), gravity causes an increase (~10%) in the 
heat transfer on an elliptical tube compared toa circular tube of the same surface area 
due to an increase in "effective" gravity. Under conditions of forced convection, both 
shear stress approximations indicate a small decrease (<< 2%) in performance due to the 
reduction in interfacial shear as a result of the better streamlined shape of the elliptical 
tube when compared tothe circular tube in the same free stream vapor flow. However, 
when pressure drop effects are also considered, the higher allowable vapor velocity over 
an elliptical tube (for the same pressure drop as a circular tube) results in an increase 
in the mean heat-transfer coefficient of 15-20%. For conditions of mixed convection, 
the condensate film thickness is controlled by both gravity and vapor shear effects. The 


cross-over between heat transfer improvement (gravity dominant flow) and _ heat 
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transfer reduction (vapor shear dominant flow) for the elliptical tube is dependent on 
the vapor boundary-layer separation point. For large condensation rates (high vapor 
suction), vapor boundary-layer separation is delayed and the cross-over point shifts to 
higher velocities (smaller F). For low condensation rates (low vapor suction), movement 
of the vapor boundary-layer separation point is minimal, resulting in a minor effect on 
heat transfer and a cross-over point that shifts toward lower velocities (large F). 
Pressure gradient and surface tension each lead to a small decrease in the mean heat- 
transfer coefficient (< 2%) for an elliptical tube. 

In general, therefore, the outside heat-transfer coefficient of a condenser tube 1s 
improved by using an elliptical tube geometry. Approximation of the single-phase 
inside heat-transfer coefficient using a hydraulic diameter and a Seider-Tate type 
correlation indicates that the elliptical tube also has better inside heat transfer 
performance compared to an equivalent surface area circular tube. Thus, the overall 


heat-transfer coefficient should be improved. 


85 


VI. RECOMMENDATIONS 


This study represents an initial step in determining if a shell-and-tube condenser 
using horizontal elliptical tubes has better heat transfer than a comparable circular tube 
bundle. The effects of multi-tube arrangements (i.e., the influence of adjacent tubes 
and the spacing between tubes) on vapor flow characteristics needs to be evaluated to 
determine if further improvements in the overall heat-transfer coefficient can be 
achieved. It is suspected that the interfacial shear of the condensate film on elliptical 
tubes will be further reduced but the corresponding reduction in vapor pressure drop 
across the condenser will result in a net gain in heat transfer compared to a circular 
tube bundle. An analysis similar to that conducted by Aoune and Burnside [46, 47] can 
be used to evaluate multi-tube arrangements. 

Further studies should be conducted ona single elliptical tube to better determine 
the effect of surface tension on heat transfer. The surface tension model used in this 
study was based only on the curvature of the tube surface and does not take into 
account the relatively thicker film over the lower half of the tube. The curvature of 
the actual condensate film surface should be analyzed in a manner similar to Krupiczka 
[29]. This analysis, together with the asymptotic shear stress assumption, would involve 
the solution of a second order ODE. 

Analysis of the effect of offsetting the elliptical tube at an angle, «, with respect 
to vapor flow and gravity must be conducted to determine its effect on the heat- 
transfer performance of the tube. This is important since in real condensers, vapor will 
approach the tube at various angles. Additionally, since the elliptical tube is a 


streamlined body, vapor flow in a direction other than the direction of the major axis 
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may result in a lifting force on the tube which may result in tube vibration and wear. 
These effects would have to be studied and methods to minimize any resultant vibration 
must found. The problems associated with joining elliptical tubes to the condenser tube 
sheet and baffle plates also needs to be addressed. 

Finally, experiments must be conducted on elliptical tubes to validate the 


predictions of these theoretical models. 
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APPENDIX A - DERIVATION OF TANGENTIAL VAPOR VELOCITY OVER AN 


ELLIPTICAL TUBE FROM POTENTIAL FLOW THEORY 


The basic technique used to determine fluid flow velocity over a body is to 
employ a conformal transformation to translate the body from one complex plane into 
another complex plane. From this, the velocity of the transformed body can be readily 
computed using potential theory. For this particular application, the fluid velocity 
along the surface of a circular cylinder is determined and then transformed into the 
complex plane in which the cylinder has an elliptical shape. 

Consider acircular cylinder in the complex plane, z,,as shown in Figure A-1. The 


surface of the circle is defined by 


Z, = exer (a) 
(AS) 


x+y =c?  (b) 


The conformal transform totranslate the circular body from the z, plane to an elliptical 


body in the z, plane is given by 





BeZ- = , (A-2) 
q 
thus 

2 =, +o) = a = ee (A-3) 

x, 4 ly, c? ee 

The real and imaginary coordinates in the z, plane are given by 
f : E } (A-4) 

x5 = x,|1 >= —| “Gee y, = yee eee 
ee of 
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Figure A-1. Circular Cylinder in z, Plane, Elliptical Cylinder in z, Plane. 


Substituting Equation (A-4)into Equation (A-1b) results in 


2 Z 
xy, 
= 2S ae ; (A-5) 


which is the expression for an ellipse in cartesian coordinates where the semi-major 
axis, @, 1S given by c + (1/c) and the semi-minor axis. b, by c - (1/c). The complex 


potential function in the z, plane for the circle is given by 


“i> Uc? a A-6 
w(Z,) = UL z,e a = + iv, ( ) 
2) 


where ® is the velocity field potential function and wis the stream function. The fluid 


velocity in the z, plane is determined from the complex derivative, 


> Ul cre? : 
ey gE Oy AD) 
az, ra 
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where u, and v, are the fluid velocity components in cartesian coordinates. The fluid 


velocity in the z, plane (for the ellipse) is found by the chain rule, thus, 


dw _ dw (A-8) 


Combining Equations (A-1la),(A-2)and (A-7) results in 


dw -U., c? (1 @8 - 72) _ gini2) (A-9} 
a, c2 ei2® 4 1 
The magnitude of the fluid velocity is found by multiplying the complex velocity by 


its conjugate. Using Euler’s formula, 


e® = cos8 + i sind , (A-10) 


and with some algebraic manipulation, the potential velocity of the fluid at the surface 


of the ellipse is given by 


y2 . dv dw _ U2 c* (2 + 2cos26) (A-11) 


= dz,dz, c¢4 +1 + 2c* cos20 


where the bar denotes the complex conjugate. Equation (A-11) is further simplified 


using the trigonometric identities 


cos260 = cos*6 - sin*0 


(A -12) 
cos26 = 2cos’6 - 1 
and the relationship between c, a and bD, 
4c? = (a+b) (A -13) 
Thus the potential velocity is given by 
Ce __ (+k) cosO (A-14) 


" cost + k?sin?6 


The angle, 6,is measured from the x, (horizontal) axis. The relationship between @ and 
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, the angle measured from the y, (vertical) axis, is given by 


o-—-¢ | (A-15) 


Note that is an angle defined in the z, plane for the circle and is the parametric angle 


defined in Section II.A . The potential velocity is thus given by 


u, = vu, —G+b sind sy, _(1+8 


F ysin? + k*cos* : V1 + k’cot*o 


(A - 16) 
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APPENDIX B - COMPUTER CODE FOR ASYMPTOTIC SHEAR STRESS MODE 
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x 
x 
X 
* 
X 
X 
X 
X 
x 
X 
x 
x 
x 
X 
X 
X 
x 
x 
x 
x 
x 
x 
x 
% 
% 
x 
x 
% 
% 
% 
% 
X 
% 
x 
x 
% 
x 
% 
X 
x 
x 
% 
X 
X 
X 
x 
X 
X 
X 
x 


THIS PROGRAM NUMERICALLY SOLVES THE PROBLEM OF MIXED 
CONVECTION CONDENSATION WITH SURFACE TENSION ON AN ELLIPTICAL 
CYLINDER. THE PROBLEM FORMULATION USES THE ASYMPTOTIC SHEAR 
STRESS APPROXIMATION ANALYSIS METHOD. 


THE PROGRAM USES AN AUTHOR DEVELOPED FORWARD STEPPING 
ALGORITHM BASED ON A GTH ORDER ADAMS PREDICTOR-CORRECTOR 
METHOD TO SOLVE A 1ST ORDER ORDINARY DIFFERENTIAL EQUATION. 


THIS PROGRAM CAN BE RUN WITH DIMENSIONLESS PARAMETERS 
(F,G,P, BO) SPECIFIED OR WHERE THE PARAMETERS ARE CALCULATED 
BASED ON SPECIFIED FLUID PROPERTIES. (MODIFICATIONS TO THE 
PROGRAM CREMOVAL OF UNNECCESSARY STEPS BY "COMMENTING®™ THE 
STATEMENT OUT OF THE PROGRAM) IS REQUIRED. SPEGTELTED FLUIDS 
PROPERTIES INCLUDE T SAT, T WALL, AND UINF. 


NUMERICAL RESULTS ARE DIRECTED TO AN EXTERNAL DATA FILE. 
MAJOR VARIABLES USED ARE: 


THETA PARAMETRIC ANGLE 

x STREAMWISE LENGTH 

DEFF EFFECTIVE DIAMETER 

DEL CONDENSATE FILM THICKNESS 


DELND,2Z DIMENSIONLESS FILM THICKNESS 
els 5 Jee DIMENSIONLESS PARAMETERS 


JARPRL 

UINF VAPOR FREE STREAM VELOCITY 
LS VAPOR SATURATION TEMPERATURE 
TW WALL SURFACE TEMPERATURE 
DENSL CONDENSATE DENSITY 

DENSV VAPOR DENSITY 

CONDL CONDENSATE THERMAL CONDUCTIVITY 
VISL CONDENSATE VISCOSITY 

HFG LATENT HEAT OF VAPORIZATION 
NU NUSSELT NUMBER 

NUAVG MEAN NU NUMBER 

ALPHA HEAT TRANSFER COEFFICIENT 
QFLUX HEAT FLUX 

A SEMI-MAJOR AXIS 

K ECCENTRICITY 

XFN PARAMETRIC RADIUS 

G GRAVITY FUNCTION 

Fl VELOCITY FUNCTION 

Ela DERIVATIVE, VELOCITY FUNCTION 
EZ PRESSURE GRADIENT FUNCTION 
ES SURFACE TENSION FUNCTION 


BE HE HEHE HE HEHE IEE IE IEE HE HE HE HE HE HE HE HEE HE HE IE IE HE HE HE HE HE EE HE HE HE DE DE EE HE HE FE HE DE HE EE HE HEE HEE HEE EK EH EE EK ERR KR 


PARAMETERCPI=3.141592654,GRAV=9.807) 
REAL Z2(¢0:180000),THETAC0:180000),X(0:180000),H,DEFF,X1,X2,X3, 


+NUAVG,NU,NUORE,UINF,TS,TW,TL, DENSL,VISL,TG,CONDL,HFG, B2, DELND, 
+DEL,QFLUX,A,K,RE, SUM, XFN,ALPHA,ALPHAQ, DQFLUX,F,PHI, 
+A1,A2,A3,A4,A5,A6,A7,A8,A9,A10,A11,TF,PSAT.TB,TC,TD,R,TE, DENSV, 
+SIGML,FCN,P,C1,€2,C3,rC0,FCl, Fezp FCs 7m clc,2P3,2F9, 50, 
+€4,C5,F21,F31,NUCOEFF 


INTEGER N,J,JJ,IC 
COMMON /GEOM/A,K, DEFF/PARAM/F, BO, P 
EXTERNAL FCN,XFN,G,G1,F1,F11,F2,F3 


OPEN(20,FILE="/MIXED1 OUTPUT Al',STATUS='OLD' ) 
OPENC350,FILE="/DEL2 DATA Al',STATUS='OLD') 
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Plies: 


KRR 


RRR 


RRR 
RX 


10 


RRR 


MIXED1 FORTRAN A 


INPUT FLUID PARAMETERS %*X 


PRINTX, ELLIPSE OF THE FORM (X/B)¥¥2+CY/A)**2=1' 
PRINT, *INPUT "A™ DIMENSION (M)7' 


READX%,A 

PRINT®, INPUT T SAT AND T WALL (K)7'° 
READX, TS, TW 

Pru %, Ittmet ANGLE STEP SIZE (DEGREES)?°* 
READX,H 


Paes, ‘Shee sIZE, H = ',H 
PRINT, ‘INPUT UINF?' 
READ, UINF 

PRINT,  INRUT aK? ° 

READX,K 

PRINT, ‘INPUT F7'* 

READX, F 


N=180/H+1 
H=PI/180.*H 
BO=0. 

P=0. 


DETERMINE FLUID PROPERTIES HHH HH HH 


A1=15.49217901 
A2=~-5.678371/7695 
A3=1.4597584637 
AG=13.877000608 
A5=~-80.887673591 
A6=123. 56883468 
A7=~-188 .321212064 
A&8=660.91763485 
A9=~-1382.4740091 
A10=1300.10901384¢ 
A11=-449.39571976 
TF=TS/1000. 
PSAT=1.E+6 XEXPCALFA2/TF+ASXLOGCTFI+T FRAG+ASRT FRXZ+A6 XT FRRS+A7® 


+T FXXG+A8 XT FHX5S+ADRTFRX6+AL0OXTFHX7+A11 TF 8 ) 


TB=1500./TS 

N@=275*L OG =EXP(=-TB)) 

eee ee SL BSTC SSRR EE) TB). 0005882xT5 
=461. 

TE=2.¥PSAT/CRETS) 

DENSV=TE/(C1.+SQRTC1.+2. XTD¥XTE)) 

TL=02.*TNt+TS)73. 

Dest =1.7C0012674-TL (2 . 029LSE-6-TL¥S .83353E-9)) 

VISL=H2.414E-5%10%*%0247 .8/0TL-140.)) 

EG=1TL/7275.15 
CONDL=-.92407+TGXC2.8395-TG*C1.8007-TG¥C.52577-TGX. 073449) )) 
HFG=3468920 .-TS*(5707 .4-TS¥*C11.5562-TS*.0133103)) 

SUG «- UUU oe C I baci ao, See e— 1 5o XC —27 3+ 75.6)71000. 


DETERMINE THETACI),X(CI),DEFF HH HHH HH 
USE SIMPSON'S RULE TO NUMERICALLY INTEGRATE ALONG SURFACE Xxx 


THETACO)=0. 
XC0)=0. 


DO 10, I=1,N 
THETACI)=THETACI~-1)+H 
XL=A¥XFNCTHETACI~-1)) 
X2=A¥XFNCCTHETACI~1)+THETACI))®.5) 
X3=AXXFNCTHETACTI)) 
XCIJV=EXCI-1LI+H/6 . ¥CX14+4. *X2+X3) 
CONTINUE 


DEFF=2.¥*X(N)/PI 

DETERMINE DIMENSIONLESS PARAMETERS AS REQUIRED RRR 
F=VISL*¥DEFFXHFGXGRAV/ CUINFXX2*CONDL XC TS-TW)) 
UINF=SQRTCVISLXHFGXGRAVZ (C FXRCONDL¥CTS-TW) )) 
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PILE: MIXED! FORTRAN A 


RRR 


RRR 
RRR 


50 


60 


RRR 


P=DENSVXHFGXVISL/(CDENSL *CONDLXCTS-TW) >) 
BO=SIGML/ ( DENSLXGRAVXDEFF 2) 
RE=DENSL*¥DEFFXUINF/VISL 

PRINTX, ‘COMPLETED THETA,X,DEFF' 


WRITEC20, 9900) ‘SAT. TEMP . CK) = ',TS,'WALL TEMP. = ',THW, 


+'I7bOn= "BUR P= oP ee ',F,*UINF CM/S) = °,UINF,'REL = ',RE 


WRITEC20, 5010) imam LENGTH CM) = ',A,'"K" ECCENTRICITY = ',«, 


+' DEFE Gap e-2, DEFF 


WRITEC20, 5020) 
DETERMINE INITIAL CONDITION, ZO KK KKK KK 


B2=SQRTCUINFXDENSL/CDEFFXVISL) ) 
THETACO)=1.E-9 

X(0)=1.E-9 

CALL DRVCE2Z,0.,F21) 

CALL ORVGES> 0. E31) 
Cl=(2*%FXG1(0.)+DEFF/SAXPXF21+3%( DEFF/A) XX2XFXBOXF31)/3. 
C2=F11(0.) 

C3=-2.*AXXFNCO.)/DEFF 
ZCOd=SQRTCCSQRT CC2*X2-G*C1 C35 3 -C23702.%C1)) 
ALPHAO=CONDL*B2/Z(0) 

SUM=0. 


DETERMINE DIMENSIONLESS FILM THICKNESS, Z HR KX 
ADAMS-MOULTON PREDICTOR/CORRECTOR ALGORITHM RRR 


FCO=FCNC(THETAC0),2(0)) 

ZP1=Z(0)+HXFCO 

ZCLI=ZC0)+H/2 . XC FCO+FCNCTHETAC1),ZP1)) 
FCL=FCNCTHETAC1),201)) 

ZP2=ZC 1) +H72. 405. Aer Guo 
ZC2)=ZCLI+H/12. XC S*FCNCTHETAC2), ZP2)+8 . XFC1-FCO>D 
FC2=FCNCTHETAC2),2¢2)) 

ZP3=Z( 2) 4H 2 XC 25. 4hCc 6 ener. xipeuD 
Z03)=Z020H729., C9 APENCTHEWAC@sD > Zr OD TLD Ap eZao gel tFCo) 
FC3=FCNCTHETAC 33,2033) 


DO 50,1=4,N 
ZP4=ZCI <LI 4H/29.4C55.. XFC5-59 . XECZtS/ . XRG =o cu) 
ZCI DSZCI-1L)+H7720.*C251 . XFCNCTHETACI) , 2P4)+646. ¥FC5-264.% 
FC2+106.. 4FCl=19" xEco> 
FCOSher 
FCL=FC2 
FC2=FC3S 
FC3=FCNCTHETACI},Z¢1)) 
C4=(2.*FXGCTHETACI) )+DEFF/AXPXF2¢( THETACI) )+3.*( DEFF/A) XX2*F 
XBOXFSCTHETACI) ))¥ZC1I) X¥34+FICTHETACID D¥2Z¢0T1) 
C5=C4+FICTHETACI 3) *ZC1) 
IFCC4.LE.0 3 He 
PRINTX,'CG LIMIT' 
Ic=I 
GO TO 60 
ENDIF 
IFCCS.LECO.) THEN 
PRINTX, * CSeEiMin* 
Ic=I 
GO TO 60 
ENDIF 
CONTINUE 
IC=N-1 


PRINTX, "COMPLETED Z DETERMINATION! 

DETERMINE HEAT TRANSFER PROPERTIES PES SS ES SS 
J=0 

JJ=0 


DO 100,I=1,1IC 
DELND=Z(1) 
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PREE: MLGEDL FORTRAN A 


DEL=DELND/B2 
ALPHA=CONDL/DEL 
SUM=SUM+CALPHAOXA®XFNCTHETACI-1) )+ALPHAXAXXFNC(THETACI) 2 ¥H¥.5 
ALPHAO=ALPHA 
IFCJJ.EQ.J) THEN 
QFLUX=ALPHAX¥(TS-TW) 
NU=DEFF/ DEL 
DQFLUX=1./DELND 
WRITEC20,5030) THETACI)*180./PI,XCI)*¥2./(DEFFXPI), 
ap DEL¥1000.,QFLUX/1000.,DELND, DQFLUX,NU 
WRITEC30,%) XCI)¥2./(DEFFXPI),DELND 
J=J+1.7180.%N+1 
ENDIF 
JJ=I+1 
100 CONTINUE 


WRITEC20,5030) THETACIC)*180./PI,XCIC)*2./( DEFFXPI), DELX1000., 
+QFLUX/1000.,DELND, DQFLUX, DEFF/ DEL 
WRITEC30,%) XCIC)*2./(DEFFPI)D, DELND 


NUAVG=SUM¥2/(PIXCONDL ) 
WRITEC20,5040) "NU/SQRTCRE) = ',NUAVG/SQRTCREJ,'NUAVG = ',NUAVG 


5000 FORMATC1X,T5,700'*') 77 1X,T5,*HEAT TRANSFER PROPERTIES 7/7 2¢1X, 
+T10,A17,F6.2 7 3,4€1X,T10,A17,E9.3 7)) 

5010 FORMATC 7 1X,7T5,*ELLIPSE GEOMETRY' 77 3(€1X,T10,A19,F6.4 7 J) 

5020 FORMATC 7 1X,T7,'THETA',T15,'DIM X',7T25,*DEL',T32,'HT FLUX', 
Sooo Olnmoele toe, D1M',158,'NUSS #' 7 T7,'CDEG)',1T25,'C(MM)',T32, 
OA et Outil Ue T6m7 C'=")9,715,60°S'),725,70°5')9,732, 
pee = oe = '),750,/0'="),758,6C'=*) 4) 

5030 ee i Cae once, F7 4, 150,77 .4, 
fe > : 

5040 FORMATC 7/7 3€1X,75,A17,F8.4 7)) 

6000 FORMAT(1X,F4¢.1,F10.3,F7.3) 


END 
HEHE HEHE HE HE HE HE EE DE HE HE 3E HE HE HE HE HE DE HE HE DE DE HE HE HE DE HE HE DE HE HE HE DE DE HE HE DE HE IE HE HE DE IE 3 HE DE HE IE DE DE HE HE HE HE DE HE HE DE HE HE DE HE DE HE HE HEE 
FUNCTION XFNCPHI) 
HE 3 HEHE HE HE HE HE 3E 3E HE HE HE 3 HE 3E HE HE IE 3E HE DE DE HE BE HE DE IE DE DE HE HE IE DE HE IE DE DE 3E DE DE DE IE HE DE HE HE DE DE HE DE DE HE DE HE 3E 3E HEE 3E HE DE DE HE DE HE HE 3 HE 
? THIS FUNCTION SUBPROGRAM IS USED TO DETERMINE STREAMWISE LENGTH 
* AS A FUNCTION OF THETA. XFN REPRESENTS PARAMETRIC RADIUS. 


% 
HEHEHE HE HE HE HEHE HEHE HE HE HE HEHE HE HE HE HE 3E 3 IE 3 3E HE HE HE DE HE HE HE DE DE IE DE 3K HE DE DE HE HE DE DE HE DE IE IE DE 3E HE HE DE HE HE HE HE IE DE IE HE 3E HE DE HE HE DE HE HE HEHE 
REAL XFN,PHI,A,K,DEFF 
COMMON /GEOM/A,K, DEFF 
XFN=SQRT C CKHH2+. 2580 (CKHR2-1 . HSINC2¥PHI) 12970 CCOSCPHI ) ) 2 
++(CKXSINCPHI) ) ¥¥2) ) 
END 
HEHEHE HEHE HE HE HE HE HEHE HE HE HE DE HE HE 3E HE HE HE DE DE 3E HE HE HE DE DE IE IE DE 3E HE DE DE IE HE DE HE IE DE 3E HE HE DE IE DE HE DE HE DE HE IE DE HE DE 3E IE DE 3E HE 3E HE IE 3E 3E HE 3E 
FUNCTION GCPHT) 
JE 3 IE 3 HE HE 3 HE HE HE HE HE DE HE HEE HE HE DE HE HE 3 3 HE HE DE HE HE HE HE DE HE HE HE HE HE HE HE HE HE HE DE HE HE HE HE HE HE HE HE HE HE DE HK HE DE NE HE DE HE HE 3E HE DE HE HE 3E HE 2K HE 
% 
¥ THIS FUNCTION SUBPROGRAM REPRESENTS THE GRAVITY FUNCTION. 
x 
HEH HEHEHE HEHEHE HEHE HEH HEHEHE HEE HE HEHE HE HEHEHE HE HE HE HE HE HE HE HE HEHEHE HE HE HEH HE HE HE HE HEHE HE HE HE HE HE HEE HE HE HK HE HK HE HE HE HE HE HE HEH HE HEHEHE 
REAL G,PHI,A,K, DEFF 


COMMON /GEOM/A,K, DEFF 
aera ee TNCRMY I ARARCKE COS CPHT) 18X29 


2G HEHE HEHE IE HE HE HE HE HE IE FE IE DE FE DE DE FE IE DE DE FE FE FE IE DE HE HE IE DE HE DE 3E IE HE IE IE DE IE DE DE DE 3E IEE HE DE FE DE DE IE DE DE HE HE HE HE HE HE HE HEHE HE HEH HE HHH 


43) 


FILE: MeED) FORTRAN A 


FUNCTION GICPHI) 


KR HK KHK HK KK HH KKH KKK HHH HK HK KKK KK HH KH HK HHH HHH HK HK HH EH HE EK EK KKK RRR 
< 


x THIS FUNCTION SUBPROGRAM REPRESENTS THE STREAMWISE DERIVATIVE 
xX OF THE GRAVITY FUNCTION. 
* 


PE SSE PSLESS ESS S SSS ESSE SSS ES eS ee esse see eee eee eee sees ees ese sete sess essere: 
REAL G1,PHI,A,K, DEFF 

COMMON/GEOM/A,K, DEFF 

COD ee on eae ace 


HE HEHE HE HEHEHE HE HE HEHEHE HEHE HE HEHE HE HE HEHE HE HE HE DE HE HE HE HEH DEH HEHE HE HE DE HE HE HE DE HE HE HE HE HE HE HE HE HEHE HE HEHE HE HE HEHE HE HEH HE HE HEH KEK 
FUNCTION F1CPHI) 


KKK HK HHH KKH KKK KKK KK KK HK KKK KKK KKK KKK KKH KK KKK HHH KK HK KKK KKK KK KKK KK RK RK KKKK 
% 


x THIS FUNCTION SUBPROGRAM REPRESENTS THE VAPOR VELOCITY 
x FUNCTION 
x* 


HE HK HEHE HEHEHE HE HE HEHEHE HK HEH HK HK HHH HH HK HK HHH HK HK HK HK HHH KK HHH HK HK HK HK KK KH HK RK KKK KKK KK 


REAL F1,PHI,A,K,DEFF 

COMMON /“GEOM/A,K, DEFF 

Bn ee Fe SS a Ca ce 
ND 


HEHE HEHE HEHEHE HE HE HEHEHE HEHEHE HE HE HEHE HEH HE HEHE HE HE IE HE HE HE HE DE HE HE HE HE HE HE HE DE DE DE NE HK HE HE HE HE HE HE HE DE HE DE HE DE HE HEE HE HE HEE HEE EK EEN 
FUNCTION FIICPHID 


HEH HEH HE HK HE HE HE HEHE HK HK HK HE HK HE HEHE HEHE HE HK HE HHH HEHE HK HK HE HEH HE HK HH HEH HHH MH HHH HK HHH HK KKH KK HK HK KKKKK 
X 


x THIS FUNCTION SUBPROGRAM REPRESENTS THE STREAMWISE DERIVATIVE 
x OF THE VAPOR VELOCITY FUNCTION. 
% 


HEHE HE HEHE HE HEE HE HE EK HEE HEE HE EEE KE HE HK HEE HK HK HK HHH HE HEH HE HK HK KH KKK HH KKK HHH HK KKK KKK KKK KKK 


REAL F1ll,PHI,A,K,DEFF 

COMMON/GEOM/A,K, DEFF 

FILS CL . HRIRAACECOS CPI) ZS QR IRCT SOL. 12 ah amamaiiseciecas) 1 aie 
ND 


HEH HEE HEHEHE HE HEHE HE DE HEHE HE HEHE HE HE HEHE HEE DEH HE DE HE HE DE HE HE HEHE HE DE DE HE HE DE HE HEH HE HE DE HE HE DE DENK DE HE HE HE HEE DEE HEHE HEE HE HE HEE EK 
FUNCTION F2CPHI) 


eS SSS SESE SESS TESTS TS TSS E SE SESS SSS ESS S SSE SESE SSS SSS ESE STS EST TS SST TSS SS 
R 


ES THIS FUNCTION SUBPROGRAM REPRESENTS THE PRESSURE GRADIENT 
x FUNCTION. 
X 


HEH HE HEHE HEHE HEHEHE HEH HEHEHE HEHEHE HE HEH HEH HEHE HK HH HEH HH HH HK HHH HH HEH HK HHH HK HHH KKK HK HHH KKK KKK KKK KKK 


REAL F2,PHI,A,K, DEFF 

EXTERNAL XFN 

COMMON’ GEOM/A,K, DEFF 

F2=CC1L.+K) XK) XX2HSINC2RXPHID/SCCSINCPHI ) ) X¥2+CKXCOS CPHI ) ) XX2) 
“a 


3 HE HE DE HE HE DE HE DE DE DE DE DE DE HE DE DE HE HE DE DE DEKE HE HE HE HEE HE HEE HE HE DE DE HE DE DE HEE DE HE HE DE HEHE HEHE NE HE HEE HEHE HE HEHE HE HEH HEHE HE HEH HE 
FUNCTION F3CPHI} 

HEHEHE HE EDIE IH HEHEHE HEN HEH HEHE HEHEHE EIEN IEEE IIE IEEE EE 

: THIS FUNCTION SUBPROGRAM REPRESENTS THE SURFACE FUNCTION. 


me 
HEHE HEHEHE HEHE HE HEE HE HE HE IE HE HE HE HE HE IE HE HE HE HE IE HE HE HE IE HE IE EE EE HE EE EEE EE IE KE EEK HEHEHE HK HK HK HK HK HK KKK KK 
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PEs eee D1 FORTRAN A 


REAL F3,PHI,A,K, DEFF 

EXTERNAL XFN 

COMMON/GEOMZA,K, DEFF 

FS=KRC] . -KAR2)KSINGC2HPHID/SQRT CCCSINCPHI ) ) X¥2+¢(K¥*COSCPHI ) ) ¥¥2) 
i Cllalalle 

N 


PESSSESEL SSS SSS LE SES ESES ESS SSS SSS SESS SSE SSES SSE SS LESS SSS SESS SS SESE SS SS: 
FUNCTION FCNCPHI,2Z) 

HEHE HE HE HEHE HE HE IE HE HE HE HE IE HE HE DE HE IE HE HE HE IE HE HE IE DE HE HE HE HE HE HE HE HE IE HE HE HEHE HE HE IE HE HE HE HEHE HE HE HE HE HEHE HE HE HE HE HEHE HE HEHE HE HEHE HEE EK 

. THIS FUNCTION SUBPROGRAM IS USED IN A PROBLEM-SOLVER 

x SUBROUTINE TO SOLVE THE O.D.E. 

x DZ/DX = FCPHI,2Z) 

39099 9 III III III III III 1 IRA RII III III 
REAL F2l,F31,PHI,FCN1,FCN2,FCN,A,K, DEFF,F,BO,P 
COMMON /“GEOM/A,K, DEFF/PARAM/F,BO,P 
EXTERNAL XFN,G,G1,F1,F11,F2,F3 


CALL DRVCF2,PHI,F21) 
CALL DRVCF3,PHI,F31) 


FONL=C2.*XFXGCPHI)+DEFF/SAXPXF2CPHI)+35.*CDEFF/A) XX2XFXBOXF3CPHI) ) 
XCAR STP LUC PHE 1x2 

FON2=2 . XAXXFNCPHI)/DEFF~-(2.¥*FXG1 (PHI )+DEFF/AXPXF21+3. XC DEFF/A) X¥2 
+XFXBOXF 31) XZ2X¥4/3.-FL LC PHI) XZ*X2 

FCN=FCN2/FCN1 

END 


PE SESESESESESESE SSS ES SSE SES ESL ESE SESE SSS ESS SSS SS SESE SS SL ESESESSS SSIS ESS 
SUBROUTINE DRVCFN,PHI,FNL) 

HEHE HEHE HEHEHE HE HEHE HEHE HE HE HE HE HE HE HE HE HE HEHE IE HE HE IE HE HE HE HE HE HE HE HE HE IE HE HE IE HE HE HE HE HE HE IE IE HE IE HE IE HE HE HE HE HE HE HE HEHE HEE HE EHH 

x THIS SUBROUTINE NUMERICALLY DIFFERENTIATES A FUNCTION. 

IIIS IH EH HI I I EEE EEE ERIE HII OIE 


REAL FN,PHI,FN1,EPS 

FNS SCEN Cranes! .E-6 )-FNCP HI =) [E=6))72 .E-6 
RETURN 

END 
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APPENDIX C -COMPUTER CODE FOR TWO-PHASE BOUNDARY LAYER 


MODEL 


% 


PROGRAM VAPBL1 


HH HEH HHH HH HH HHH HHH HH HHH HH HK HHH HHH HHH HHH HHH HH HK HHH HH HH HAHAH KKK NRK HK RRA KKRRKK 


THIS PROGRAM NUMERICALLY SOLVES THE PROBLEM OF MIXED 

CONVECTION CONDENSATION WITH SURFACE TENSION ON AN ELLIPTICAL 
CYLINDER. THE PROBLEM FORMULATION USES THE SOLUTION OF THE 
TWO-PHASE BOUNDARY LAYER EQUATIONS FOR THE CONDENSATE AND VAPOR. 


THE PROGRAM USES AN IMSL NUMERICAL PROBLEM SOLVER, DIVPAG, 
TO SOLVE A 1ST ORDER ORDINARY DIFFERENTIAL EQUATION. THE 
PROBLEM SOLVER IS WRITTEN IN DOUBLE PRECISION WHILE THE MAIN 
PROGRAM IS IN SINGLE PRECISION. THE MAIN PROGRAM MUST BE 
COMPILED IN DOUBLE PRECISION (AUTO DOUBLE). 


THIS PROGRAM CAN BE RUN WITH DIMENSIONLESS PARAMETERS 

(F,G) SPECIFIED OR WHERE THE PARAMETERS ARE CALCULATED 

BASED ON SPECIFIED FLUID PROPERTIES. MODIFICATIONS TO THE 
PROGRAM CREMOVAL OF UNNECCESSARY STEPS BY "COMMENTING" THE 
STATEMENT OUT OF THE PROGRAM) IS REQUIRED. SPECERIED FeUrD 
PROPERTIES INCLUDE T SAT, T WALL, AND UINF. 


IF THE ALGORITHM HAS DIFFICULTY LOCATING THE VAPOR SEPARATION 
POINT, THEN THE PROGRAM NEEDS TO BE MODIFIED TO INPUT THE 
PARAMETRIC AND, PHI, FOR THE SEPARATION POINT. 

NUMERICAL RESULTS ARE DIRECTED TO AN EXTERNAL DATA FILE. 


MAJOR VARIABLES USED ARE: 


PHI PARAMETRIC ANGLE 

X STREAMWISE LENGTH 

DEFF EFFECTIVE DIAMETER 

DEL CONDENSATE FILM THICKNESS 
DELND DIMENSIONLESS FILM THICKNESS 
é DIMENSIONLESS VAPOR MOMENTUM THICKNESS 
F,JARPRL DIMENSIONLESS PARAMETERS 

UINF VAPOR FREE STREAM VELOCITY 

ibs VAPOR SATURATION TEMPERATURE 
TW WALL SURFACE TEMPERATURE 
DENSL CONDENSATE DENSITY 

DENSV VAPOR DENSITY 

CONDL CONDENSATE THERMAL CONDUCTIVITY 
Vise CONDENSATE VISCOSITY 

VISV VISCOSITY OF VAPOR 

HFG LATENT HEAT OF VAPORIZATION 
NU NUSSELT NUMBER 

REV VAPOR REYNOLDS NUMBER 

RE TWO-PHASE REYNOLDS NUMBER 
NUAVG MEAN NU NUMBER 

ALPHA HEAT TRANSFER COEFFICIENT 
QFLUX HEAT FLUX 

A SEMI“-MAJOR AXIS 

K ECCENTRICITY 

XFN PARAMETRIC RADIUS 

G GRAVITY FUNCTION 

Fl VELOCITY FUNCTION 

Fad DERIVATIVE, VELOCITY FUNCTION 
Ez PRESSURE GRADIENT FUNCTION 

FS SURFACE TENSION FUNCTION 

TAU INTERFACIAL SHEAR 

KAPPA PRESSURE GRADIENT FORM PARAMETER 
KAPPAA SUCTION FORM PARAMETER 


HEH HH HHH HK HHH HHH HEHE HEH HEHEHE HEHEHE HEHEHE HEHE HEH I HE HEHE HEHE HEHE HE HE I HEH HE HH HHH HH HHH HH HK HK HH RNR AK RK KX 


PARAMETERCPI=3.141592654,GRAV=9 .807 ,NEQ=2, NEQ1=1,NPARAM=50) 
REAL PHIC0:180000),DEL(0:180000),X€0:180000),A,TS,7W, 


+H, Al1,A2Z,A5,A4,A5,A6,A7,A8,A9,A10,A11,TF,PSAT,TB,TC,TD,R,TE, DENSV, 
+TL,DENSL,VISL,VISV,TG,CONDL,HFG,SIGML,F,K,X1,X%2,X5, DEFF,UINF, 
+REV,RE,JARPRL, DELOO, DELO,Z0, TAU, DTAU, FDELO, DFDELO, DEL1,21,ALPHA0, 
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Pebe:  VAPBL FORTRAN A 


HR 


KR 
RRR 


+SUM,AAC1,1),FCN,FCNJ,HINIT, PARAMCNPARAM) , KAPPA, KAPPAA, DKAPPA, 
+TOL,XEND,XBEG,YC2),DDEL,DZ,Y1€1),NU,NUAVG,NUORE, DKAPPAA, XO, Dx0, 
+DELND, DELL,ALPHA, QFLUX, DQFLUX, D2UPHI, DUPHI, UPHI 

INTEGER N,JJ,II,IDO,IMETH 

COMMON /GEOM/A,K, DEFF/PARAM/ JARPRL,F 


Sa ZFN,DELFN,XFN,G,G1,F1,F11,NUSSELT, DRV,FCN,FCNJ,IVPAG, 
+33 


OPENC20,FILE="/VAPORI OUTPUT Al',STATUS='OLD') 
OPENC30,FILE="/VAPOR2 OUTPUT A1',STATUS='OLD') 


PRINT, *ELLIPSE OF THE FORM (X/B)¥¥2+(CY/A) ¥X2=1' 
PRINTX, ‘INPUT "A™® DIMENSION (MM)?! 


READX,A 
PRINTX,* INPUT T SAT AND T WALL (KD)?! 
READ®, TS, TW 
ae ala ANGLE STEP SIZE (DEGREES)?! 
ADX® 
PRan Aveeomer SIZE, H = *°,H 
PRINTX,*INPUT PHI AT SEPARATION? ' 
READX¥,PHIS 
PREG, *INPUT UINF? ° 
READX, UINF 
PRINTX, ‘INPUT K7' 
READX,K 
PRig0x, “INPUT G7° 
READX, JARPRL 
N=180/H+1 
PRINT, 'N-1L=',N-1 
H=PI/180.XH 


DETERMINE FLUID PROPERTIES HHH HHH 


A1=15.49217901 

A2=-5.67837176953 

A3=1.4597584637 

A4=13.877000608 

A5=-80.887673591 

A6=123 .56883468 

A7=~-188 .321212064% 

A8=660.91763485 

A9=-1382.4740091 

A10=1300.1040184% 

A11=-449 .39571976 

TF=TS/1000. 
PSAT=1.E+6XEXPCAL+A2/TFtAS XL OGCTF)+TFRHAG+A5 XT FRX2+A6 XT FRXS+A/® 
+T FRRG+AS8 RT FRXStAIRT FRG +ALL ORT FRR7 +AL LRT PRS D 

TB=1500./TS 

IC=25et0GC1.-EXPC=1B) ) 

em TE ARTE) TB). 00048 82KTB 
R=461.51 

TE=2.*PSAT/CRETS) 

DENSV=TE/C1.#5QRTC1.4+2.%TDXTED) 

TL=02. ¥TWNtTS)73. 

MEN Salat OOl767 S—10 402. 02915E=6-TL¥3.8333E-9') D 
VISL=2.419E-5%10%%( 247 .8/(CTL-140.)) 

VISV=—-4% .978415E-6+7TS¥C5.O216E-8-1 .579E~-11XTS) 

G= IE 729-5 15 
CONDL=~.92407+TGXC2.8395-TGRC1 .8007-TGXC .52577-TG.07344))) 
MEG=5465920~15%05707 .4-1SXC11 .5562-TSX.0133103)) 

Some =( =. OU Tiee7 3, LS Ixee— eee CT=—275.15)+75.6)71000. 


DETERMINE PHICI),XC€1I),DEFF KKK HK KKK 
USES SIMPSON'S RULE TO NUMERICALLY INTEGRATE ALONG SURFACE RX 


PHICO)=0. 
X€0)=0. 
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PIECE: 


10 


HRKK 


RRR 
RRR 
RK 
KKK 
XK 


20 
30 


VAG OE FORTRAN A 


Ne et 


+PHICI))%.5) 


6.*CX1+4. *¥X2+X3) 
CONTINUE 


DEFF=Z2. <xXCND7PL 
PRINT®, ‘COMPLETED DE DETERMINATION, DE=!,DEFF 


DIMENSIONLESS PARAMETER DETERMINATION KKK 


F=VISL¥DEFFXHFG/ CCONDL XUINFX¥2*( TS-TW) )D 
REV=UINFXDENSV*¥DEFF/VISV 

RESUINFXDENSLXDEFF/VISL 
JARPRL=CONDLXCTS-TWI/HFGXSQRT CDENSLZ CDENSVXVISLXVISV) ) 
PRINT%, 'COMPLETED PARAMETER DETERMINATION! 


WRITEC20,5000) 'SAT. TEMP. (K) = ',TS,'WALL TEMP. = ',TH, 
tf = ',F,'UINF (M/S) = ',UINE,'G = ',JARPRE 

WRITEC20,5010) '"A" LENGTH (M) = ',A,'"K™ ECCENTRICITY = ',K, 
TIDEFE eG) = "HDEEF 


DETERMINE INITIAL VALUES OF Z AND DEL KKK 
USES INTIAL VALUE OF DEL FROM SHEKRILADZE-GOMELAURI XX 
FOR MIXED CONVECTION AS STARTING POINT TO FIND Z THEN Xxx 
REDETERMINES DEL FOR THE PREDECTED VABUE URSZeUSTIHG HX 
NEWTON'S METHOD, CHECKS FOR CONVERGENCE AND ITERATES xx 


DELOO=SQRTC.75*XK/FRCSQRT CC CL. tK)/K) €X2416 . XFRXA/ C3. XDEFF)) 
+=(1. +7 


DELO=DELO0 
Z20=(-.07355*XFNCO.)XJARPRL/DEFF/F11C0. I7DELO+SQRT CC .07355% 
+XFNCO.)*¥JARPRL/DEFF/SFIICO. )/DELO)*¥2+.07355*XFNCO.)/DEFFZ 
+F1160.)9))%%2 


UPHI =FrCo> 
DUPHI=F11(0.) 

CALL DRVCF11,0.,D2UPHI) 
XO=XFNCO.) 

CALL DRVCXFN,0O.,Dx0) 


KAP PA=DEFF/X0XDUPHI ¥Z0 
KAPPAA=.0682+.174¢*JARPRL¥SQRT(Z0)/ DELO 


DKAPPA=DEFFXZOXD2UPHI/X0-DEFFXDUPHIXZOXDXO0/ X02 
DTAU=6 .44%SQRT CKAPPAXKAPPAA+KAPPAA®X2 ) ¥DUPHI/SQRT(Z0)+3.22*UPHI 
+¥KAPPAAXDKAPPA/SQRT(Z0)/7SQRT CKAPPAXKAP PAA+KAPPAAXX2 ) 


FDELO=X0/DEFF-F/35.XG1(0. )*DELOX*4-DTAUXDEL 0%%3/(04.*JARPRL ) 
DFDELO=-4. XF¥G1(0. )*¥DELO*¥€3/73.-3.*DTAUXDELOX*2/(04. *JARPRL) 
DELI=DELO-FDELO/DFDELO 

IFC CABSCDELI-~DELO).GT.( .OOLXDELO. JR AND. (FDELS Gi .1.E-7)) (iis 


+ 


DELO=DEL1 
GO TO 30 
END TE 


Z1=(-.0735%*XFNCO. XJARPRL/ ( DEFFXF1100.)*DELI)I+SQRTCC.0755% 
+XFNCO. *XJARPRL/CDEPFHFIICO. UXDELI 4X24. 07 SOXXPNCO. D7 (DERE 
+PL1.0.))2) RZ 


IF (CABSC(DELI-DEL00).GT.¢.O01¥*DEL00)).OR.CABS(Z1-Z0).GT.C.001% 
+Z0))) THEN 
20=21 
DELOO=DEL1 
DELO=DELI 
GO TO 20 
ENDIF 


DELCO) =SDEL1 
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FIRECa VAP BI FORTRAN A 


RRR 
XXX 


RRR 


100 
190 


RHR 
RRR 
RRR 


200 


KR 


ALPHAQ=1./DEL(0) 

SUM=0. 

PRINTX, ‘COMPLETED 2¢0),DELC0) DETERMINATION’, '! aE a Ge 
+'DEL=', DELCO) 


EVALUATION OF ZCI) AND DELCTI) XXXX 
USES IMSL ROUTINE DIVPAG KKK 
HINIT=1.E-4 

IMETH=2 


CALL SSETCNPARAM,0.0,PARAM,1) 
PARAMC1)=HINIT 
PARAMC1L2)=IMETH 
XBEG=PHICQ) 
YC1)=Z1 
YC2)=DEL(C0) 
Lobe... E=35 
IDO=1 
DO 100,I1=1,N-1 
PeCkwEor).Gi.PHIS) GO TO 190 
XEND=PHICTI) 
CALL DIVPAGCIDO,NEQ,FCN,FCNJ,AA,XBEG,XEND, TOL, PARAM, Y) 
DELC I) =YC2) 
KAPPA=DEFF/XFNCPHICI))XFLILCPHIC I) D¥YC1) 
KAPPAA=.0682+.174¢¥JARPRLXSQRTCYCLIIZYC2) 
TAU=6 .44%F1LCPHICI)) ¥SQRTCKAPPAXKAPPAA+KAPPAAXX2)/SQRTCYC1)) 
WRITEC30,%) "PHI=',PHICI)X180/PI, 'TAU=', TAU 


its 

CHECK FOR VAPOR BOUNDARY LAYER SEPARATION X RKX 
IFCTAU.LE.0.) GO TO 190 

CONTINUE 


PRINT, *SEPARATION ANGLE AT PHI=!,PHICII)¥180./PI 
PRINTX, 'TAU=',TAU, 'DEL=!,DELCII) 


ID0=3 
CALL DIVPAGCIDO,NEQ,FCN,FCNJ,AA,XBEG, XEND,TOL,PARAM,Y) 


USES A SIMPLE NUSSELT ANALYSIS TO COMPUTE FILM THICKNESS *¥ XXX 
DONNSTREAM OF THE SEPARATION POINT. PROBLEM SOLVER IS ¥XXX 
IMSL DIVPAG ¥XXX 


PF CI, NEWON= Oe THEN 
XBEG=PHICII) 
IDO=1 
YIC1)=DEL CIID) 
DO 200, I=II+tl,N-1l 
XEND=PHICT) 
CALL DIVPAGCIDO,NEQ1,NUSSELT,FCNJ,AA,XBEG, XEND, TOL, PARAM 


pale) 
DELCI)=Y1C1) 
CONTINUE 
ENDIF 


IDO=35 
CALL DIVPAGCIDO,NEQ1,NUSSELT,FCNJ,AA,XBEG,XEND, TOL, PARAM, YL) 


PRINTX, ‘COMPLETED ZCI), DELCI) DETERMINATION’ 
WRITEC20,5050) PHICII)*130./PI 
WRITEC20,5020) 


DETERMINE HEAT TRANSFER PROPERTIES PES SELLS SS 


JJ=0 

DO 300,1=1,N-l 
DELND=DELCI) 
DELL=DELND¥SQRTCVISLXDEFF/ CDENSLXUINF) ) 
ALPHA=1./DEL(CI) 
SUM=SUM+CALPHAOXXFNCPHICI-1L) )+ALPHAXXFNCPHICI)) )¥HX.5 
ALPHAO=ALPHA 
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FILE: SAB 8 L FORTRAN A 


IFC CJJ.E9. 10 OR-CJJ-60.0)) nen 
QFLUX=ALPHAXCTS-TW) 
NU=DEFF/DELL 
DQFLUX=1./DELND 
JJ=JJ+10./130.¥*N+1 
WRITECZO,5030) PHICIO*1 30781) G27 COELEXP Ie DELLE 
a ¥1000,QFLUX/1000, DELND, DQFLUX,NU 
ENDIF 
300 CONTINUE 


WRITEC20,5030) PHICN-1)%180./PI,X(N-1)¥2./(DEFFRPI),DELL¥1000., 
+QFLUX/1000.,DELND, DQFLUX, DEFF/DELL 
WRITEC30,%) XCN-1)¥2./(DEFFXPI), DELND 


NUAVG=SUM¥2/ CPIXDEFF) 
WRITEC20,5040) 'NU/SQRTCRE) = *,NUAVG/SQRTCRE) 


5000 FORMATCIX,T5,700'*") “7 1X,T5,*HEAT TRANSFER PROPERTIES! 77 2¢1X, 
+T1O,Al7,F6.2 / ),5C1X, 010 A Eo oe 

5010 FORMATC 7 1X,T5,'ELLIPSE GEOMETRY! 77 3¢€1X,T10,A19,F6.4 7 J) 

5020. FORMATS /_ 1X,177,'THETA', [15,80 0M xX", 125, UEC fs2, Hl FLU 
+T91,"DIM DEL',T52,'DIM',T58,"NUSS #° 7 T7,*CDEG)",1T25," CMM) "poe 
*'CKW/MZ)',TS50,°HT FLUX! 7 T6,7C€'="),715,6€"="),7235,7C' ="), tee 
#70'3'°),T41,70'='),T50,70'='),T58,60'=') 7) 

5030 SSG ye cemneegiee Tan taro eeacgae 7 - 97a 
+ ,F9. 

5040 FORMATC 7/7 3¢€1X,T5,A17,F8.4 7)) 

5050 FORMATC1X,1T6,'VAPOR SEPARATION AT PHI = ',F6.2) 


END 
PES SSSSELELE SESE S SSSESSESESS SESE SESE SSEESESES ESS SESES SES SES ESE SESE TSS ESSE SF 
FUNCTION XFNCPHI) 


SS 2.2,. 2. 2,0,.2,0,2,5, 2,2, 5,2,2,2,0,5,5,2,5,5,5,5,5,5,2, 2,2, 5,5,2,5,5,2,5,2,5,5,2,5,5,5,5,5,5,5,5,2,0,5,5,5,2,5,5,5,5,5,5,5,2,5,5, 5,5, 5,545 
¥% 


x THIS FUNCTION SUBPROGRAM IS USED TO DETERMINE STREAMWISE LENGTH 
x AS A FUNCTION OF PHI. 
x 


So 2.2.2. 2.2,2,2,0,5,.5,5,5,. 5,5, 5.2, 5,2, 5,5, 5,2,2,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,2,5,5,5,5,5,5,5,5,5,5,5,5,5,5,2,8,5,5,9,5,5,5,0,5, 545 


REAL XFN,PHI,A,K,DEFF 
COMMON /GEOM/ZA,K, DEFF 
XFN=AXSQRTC CK¥¥2+. 25%CCKR2-1. DXSINC2*PHI) ) ¥¥2)/CCCOSCPHI) )¥¥2 
eae 
ND 


3€ 3 IE 3E 3 HE DE DE BE 3 3 IEE BE 3 BE IE 3 3 HE DE BE DE DE BE BE BE DE DE BE HE DE HE HE DE HE BE DE DE 3 HE HE 3 DE DE DE DE BE BEE DE 3K 3 IE DE DE IE DE DE HE DE DE DE HE HE DEE EE 
FUNCTION GCPHI) 

DESSSESESESS SESS S SSS SESS SS SSS ESS S SESE SES ESS SE SES ESSE S ES SESE se sesesss sds +: 

GRAVITY FUNCTION 


HE HE HE HE HE HE HE DE DE HE DE IE HE DE DE DE SESE HE HE BE HEHE DE HE HE DE IE DE HE HE HE EE HE HE IE DE DE DE DEE HE DE SE DE DE DE IE HE DE HE IE HE HE HEE HE EE EEE HK HK RN 


REAL G,PHI,A,K, DEFF 
COMMON /“GEOM/A,K, DEFF 
Cae ne ne anes OSE HL J Cake el eo 
ND 
HEHE HEHEHE HE HEHE HEHE HEHE 3 HEHE DE DEI HE HE 3 HE IE HE DE HE IE IE HE DE IE HE HE DE HE 3 DE DE DE HE DE HE DE DE HE DEE 3 HE DE DE DE 3 HE DE HE DE DE DE DE HE DE HE DE HE DE 3K IE HEE 
FUNCTION G1CPHI) 


HEHE HEHE HEHE HE HE HE HE HE HE HEHE HE DED KE DE DE DE HE HE HE HE DE HE HE DE HE HE HE HE BE HE HE HE HE HE HE HE DEE HE DE FE HE HE HE SE IE HE DE EE HEE HE HEE HEHE HK HK MRK 


X 
xX STREAMWISE DERIVATIVE OF GRAVITY FUNCTION 
X 
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Powe VAP OL FORTRAN A 


KH KH KH KK HH HHH HHH HHH HH HH HK HHH HK HHH HH HHH HK HHH HHH HHH HHH HHH HHH HK HHH HK KKK KKK 


REAL Gl,PHI,A,K, DEFF 
COMMON/GEOM/A,K, DEFF 
a reeds SORTER Et RXCOSCPHL) 8X2) RRS D 
HEHE HEHE HEHE HE HE HE HEHE HE HE HE HE HE HEHE HEHE HE HE HE HE HE HE HEHE HE HE HEE HE HE HE HE HE HE HE HE HE HE HE HE HE HE HE HE IE HE HE HE HE IE HE HE HE HE DE HE HE HE HE HE HE HEHE HEHE 
FUNCTION FILCPHI) 
HEI HEHE HE HEHE HE HE HEHE HE HE HEH HE HE HE HEHE IE HE HEE DE HE HEHE HE HE DE HE HE DE HE IE DE HE DE HE HE HE HE HE I HE DE HE I DE DE I DE DE HE DE DE DE DE HE IK DE HEHE EK EN 
X 


* VAPOR VELOCITY FUNCTION 
% 


HHH HH HH HHH HHH HHH KH HK HH HHH HK KK KH KKK KKK HK KKK KK KKK KH KK HKAK KK KKRKRKK KKKKKKKK KKK 


REAL F1l,PHI,A,K, DEFF 
COMMON /GEOM/A,K, DEFF 
aia ans eR em ECOS LENT) EXC) 


HHH HHH HH HH HHH HHH HHH HH HH HHH HHH HK HH HK HH HHH HK HK HHH HH HHH KKK HH KKK HK KKK KK KAKKK 


FUNCTION FL1CPHI) 


HEH HHH HHH HH HH HEH HEH HEH HEH HEHEHE HHH HH HHH HHH HHH HH HK HH HHH KKK KKK HH HAH HK KKK AKHK KKK KKKKR 


X 
* STREAMWISE DERIVATIVE OF VELOCITY FUNCTION 
% 


HHH KH HHH HHH HH HHH HH HHH HHH HHH HH HHH HHH HH HH HK HH HHH KKH HHH KK HKHK KKK K KK KK KK KKK 


REAL F1ll,PHI,A,K, DEFF 
COMMON/GEOM/A,K, DEFF 
a erect ene eerecin COS CPI) 2X2) 885) 


HH HHH HHH HH HHH HH HK HHH HK HH HHH HK HK HHH HE HHH HK HHH HHH HHH HHH KKH KKK HHH HHH KK KHK KKK KKK 


FUNCTION DELFNCPHI,2Z, DEL, DZ) 


22.2.2. 2, 2,2, 22.022 22 222 22 22 2 2S EES PES S25 55,05, 5,5, 0,5, 5,0, 0,5,5, 0,5, 5,5, 0,5, 0,0, 5,5, 0,5, 5,5, 55,5, 5,5, 5,0,5,: 


x 
% STREAMWISE DERIVATIVE OF CONDENSATE FILM THICKNESS 
Xx 


HH HHH HH HH HH HHH HH HHH HH HHH HHH HK HH HHH HHH HK KKK KKK HK KK KK HK KK RKKK KKH KKK RKRKKKKKR 


REAL DELFN,PHI,Z,DEL,DZ,DELFNI, DELFN2, DELFNS,X,FF1,FFL1,FF12,DXFN 
+,A,K, DEFF,JARPRL,F,KAPPA, KAPPAA, XX, FKAPPA 


EXTERNAL XFN,F1,F11,G,Gl 


COMMON/GEOM/A,K, DEFF/PARAM/ JARPRL, F 


IFCPHI.EQ.0.D+0) THEN 
aaa ale 


ELS 

FFL=FICPHI) 

Pe lyker ll (PHI 

CALL DRVCF11,PHI,FF12) 

A=XFNCPHI ) 

CALL DRVCXFN,PHI, DXFN) 
KAPPA=DEFFXRFF11¥Z/X 

KAPPAA= .0682+.174*JARPRL¥SQRT(Z)/DEL 
AX=KAPPAAXKAPPA+KAPPAAXX2 
FKAPPA=SQRTCXX) 


DELFNL=.805X¥KAPPAAXDEFFXFFLXFFLILXDEL¥¥3/( JARPRLXFKAPPAXX¥SQRT (CZ) ) 


++(7 .0035E-2)%(2. XKAPPAA+KAPPA) XFFLXDEL X27 ( FKAPPAXZ)—.805*XFKAPPAX 
+FFILXDEL €X3/( JARPRL ¥ZXX1 . 5) 
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-FiLE-sVARBE FORTRAN A 


DELFN2=FXGC PHI ) X¥DEL¥¥3+3.22*FFIXFKAPPAXDEL €¥2/ ( JARPRL¥SQRT(Z) ) 
+-.14007%02.*KAPPAA+KAPPA) XFFIX¥DEL/FKAPPA 


DELFN3=X/ DEFF-FXG1 (PHI) XDEL¥¥4/3.-1.61*FKAPPAXFFIIXDEL¥%3/( JARFRL 
+¥SQRTCZ))-. 805XKAPPAAXDEFFRFFIXFFI2€SQRT (CZ) XDELX¥S/ ( JARPRLXFKAPPAX 
Ce eee a ee tA ieee a 
Xxx 

DEVEN-=@DELFNS=DELFNIXDZ)7~DELENZ 

EN DLE 

END 


HEHEHE HK HEHE HEHE HE HK EE EE HEE HEE HE HEE HEE HE EE HEE HK EE EE KE EE HE EE EE HK EE EEE EEK EEK EHH RK KK KKH KKK 


FUNCTION ZFNCPHI,Z, DEL) 


HEHE HE HE DE HE DE HE HE HE HE DE HE DE HE HE HE HE HE HE HE DE DE IE HEE DE DE DE HE DE HE HE HE HE HE DE HE DE HE DE HE DE HEHE HE HE HEHE HE HEHE HEH HK HER HEHE HK EK HK HEN 
* 

€ STREAMWISE DERIVATIVE OF VAPOR MOMENTUM THICKNESS PARAMETER 

x 


HEHE HE DE HE DE HE HE DK HE HE EH HE HE HE DK HE HE DE HE DK HE DK DE HE HEHE DE HE HE HE HE DK HE HE DEH DE HE DE DE DE HEHE DE HE DE HE HE DEE HK DE DE HE DE HE DE DE HE DEE EK EK HE 
REAL ZFN,PHI,Z,DEL,DEFF,A,K, JARPRL,F 
EXTERNAL XFN 
EXTERNAL Fl 
EXTERNAL F1l - 
COMMON/GEOM/A,K, DEFF/PARAM/ JARPRL,F 
TECPHL cQ so. Dru rine 
ZFN=0.D+0 
ELSE 
ZFN=C€ .441%C1.-2.*JARPRL¥SQRT(Z)/DEL)—-6 .XDEFFRFLICPHIDXZ/XFNCPHTI)) 
+EXENCP NL 7 Pern DEP Fr 
ENDIF 
END 


HEHE HEHE HEE HE DE HE HE IE IE IEE HE IE DE DE HE EE HE HEE HE HE HEE HE HE IE IEE IE IEE IE IEE DE IE HE IE HE HE EE IE EK EEE EEK EEE KHER KR KKK 


SUBROUTINE NUSSELTCNEQ,PHI,Y,YPRIME) 


HEHEHE HE HEH HEHE HE HEHE HE HEHE HE HEHE HE HEHE HE DE HE HEHE HE HEHE DE HE HE HE HE DK HE HK HE HE HE HE DK HEHE HEHE DE HE DE HE DE HE HEHE DE HEHE DE DE HEE DEE HE DE HE HEE EK 
x 

Bb STREAMWISE DERIVATIVE OF CONDENSATE FILM THICKNESS USING 

x A NUSSELT ANALYSIS DOWNSTREAM OF VAPOR SEPARATION POINT. 

x 

HH HH HEHE HE HE HE HE HE HEHEHE HE HEHE HEHE HE HEHEHE HE HEHE HEHE HE HEHEHE HE HEHEHE HEH HEHE HK HEHE HE HE HEHE HEHE HE HE HEHE HE DE HEHE HE HE DE HE DE HE THEE HE HEME 


INTEGER NEQ 
REAL PHI, YCNEQ),YPRIMECNEQ),A,K,DEFF,JARPRL,F 


EXTERNAL XFN,G,G1 
COMMON/GEOM/A,K, DEFF/PARAM/JARPRL, F 
YPRIMEC1)=CXFNCPHI)/DEFF-F/3. *YC1)**4*G1 (PHI) I/CFXGC PHI) XYC1) X¥3) 


RETURN 
END 


FE IE HE HE HEHEHE HE HE DE HE HE HE HE HE DK HE HE HE HE HEHE HE HE HE HE HE HK HE HE HE DE HE HEHE HEHE DE HE DE HE HK DE DK DE DE HE HE DE HE IE HE DE HE DE HE DE DE I HE DEE HK HEE EEX £ 
SUBROUTINE DRVCFN,PHI,FNI) 

HE HEHE HEHE HE HE HE DE HEE HED DE HE DE HE HE DE HE DE DE DE DE HE HE HE DE HE HE HE HE DE HE HE HE HE HE HE DE HE DE HE HE I DE DE HE DE HE 3E HE DE DE 3 DE DE HE DE DE DE HEE HE HE HE HEN 

SUBROUTINE TO NUMERICALLY DIFFERENTIATE A FUNCTION 


HE HEHEHE HE IE HE HE IE HE SE HE DE HE HE HE EE IE IE DEE DE DEE SE DE SE IE SEE DE IE DEE DEE HEE HE HE HE HE HE HE EE EK HEHE EK KE HE HK HEHE HHH HRM KKK 
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PIE: VAPBL FORTRAN A 
Reale EN, Fal, FN1,EPS 
Piieatemeritit) .—E-6)-FNOPHIT—] .E-6))/2.E-6 


RETURN 
END 


HEI HE HE FE HE HE HE HE HE HEHE HE HE 3 3 HE HE 3 DE HE 3 DE DE 3 DE DE 3E DE DE DE DE DE DE HE DE DE DE DE 3E HE 3K DE DE DE DE HE IE IE HE IE DE HE DE 3 DE DE DE DE IE 3 HE DE IE HE HE HE HE 
SUBROUTINE TAUFNCPHI,Z,DEL,DZ,DDEL,TAU,DTAU) 


HHH HH HHH HK HHH HHH HHH HH HH HH HHH HEHE HEHE HEHEHE HEHEHE HE HEHEHE HE HE HEHE HEHE HE EE HE HE HEHE HE EE HE HEE EE RHE 
x 


x SUBROUTINE TO DETERMINE THE INTERFACIAL SHEAR STRESS AND 
* ITS STREAMWISE DERIVATIVE. 
x 


BE HE HE HE HE HE HEHEHE HEH HEHE HE HEHEHE HEHEHE HEHEHE HE HE HEHE HEHEHE HEHEHE HE HE HE HE HEHEHE HHH HHH HE HEH HH HH HH HH HHH HK KH RK RR 


REAL PHI,Z,DEL,DXFN,DZ,DDEL,KAPPA, KAPPA1,KAPPAA, DKAPPA, 
+DKAPPA1, DKAPPAA, TAU, DTAU,A,K, DEFF, JARPRL,F 


EXTERNAL XFN,F1,F11 
COMMON/GEOMZA,K, DEFF/PARAMZ JARPRL, F 


CAEL SORV CELL PHI, F12) 
CALL DRVCXFN, PHI, DXFN) 


KAPPA=DEFFRXFLICPHI)*Z/XFNCPHI) 

KAPPAL=JARPRLXSQRTCZ)/DEL 

KAPPAA=.0682+.174¢%KAPPAI1 

DKAPPA=DEFF/XFNC PHI )¥(ZHFI2+F1l 1 CPHI)DXDZ)-DEFF/ CXFNCPHI) ) 2x 
+FLICPHI)XZ¥DXFN 

DKAPPAL=JARPRLXC . 5XDELXDZ-Z¥DDEL I/CSQRTCZIXDEL R¥2) 

DKAPPAA=.174%DKAPPA1 


TAU=6 .44¢%F11 CPHI)*SQRTC CKAPPAAXKAPPA+KAPPAAXX2)/Z) 

DTAU=6 .44%( ( CKAPPAAXKAPPA+KAPPAAX¥2) XFLICPHID +. 5¥*F1CPHIDX 
+CKAPPAAXDKAPPA+( 2. XKAPPAA+KAPPA ) XDKAPPAA) J/SQRT (ZC KAPPAAX 
+KAPPA+KAPPAA®#2) )-. 5X CKAPPAAXKAPPA+KAPPAAX#2) XF1CPHI) ¥DZ/ 
+SQRTCZHX SX CKAPPAAXKAPPA+KAPPAAX¥2) ) ) 


RETURN 
END 


HEE HE HEHE HEHEHE HE HE HE HE HE HE HE HE HE HE HE DE HE HE HE HE HE HE HE DE HE HE DE HE IE IE DE HE HE DE HE HE HE HE DE HE HE DE DE HEE DE DE DE HE DE DE DE DE DE DE DE HE DE DE DE DE DE EE EE 
SUBROUTINE FCNCNEQ,X,Y,YPRIME) 


HEH HEHEHE HEHEHE HEHE HE HEHE HEHE HEHEHE HEHE HEHE HEHEHE HE HEHEHE HEHEHE HE HEHEHE HEHE HH HEH HEH HHH HHH HH HHH HHH HH HH KKK HK KKK 
x 


x SUBROUTINE TO INPUT DIFFERENTIAL EQUATION TO BE SOLVED 
% By le PROBBEM SOLVER. 
x 


HEHE HE HEHE HE HE HE HE HE IE DE HE DE HE HE DE HE DE 3 HE DE HE DE DE HE DE HE DE HE DE DE DE DE DE DE HE DE HE DE DE DE DE HE DE BE DE 3 HE DE HE DE DE HE DE BE DE DE HE DE HE DE DE DE DE HE 3E IEE 
INTEGER NEQ 
REAL X,YCNEQ),YPRIMECNEQ) 
EXTERNAL ZFN,DELFN 


ebm l )=ZENGX. YC 1), YC2)) 
YPRIMECZ)=DELPFNCX,YC1),YC2),YPRIMECI)) 


RETURN 
cND 


HEHE HEHE HEHE HEE HEHEHE HE HEHE HE HE HE HEHE HE HE HEHE HE HEHE HE HE HEHE HE HE HEHEHE HEHEHE IE HE HE HE DE HE IE IE HE IE HE HE HE DE DE DE DE DE DE DK DE DE DE DE DE HE HE HEHE HE EK 
SUBROUTINE FCNJCNEQ,X,Y,DYPDY) 


HEHE HEHE HE HE IE IE HE DE DE HE DE HEHE HEE HE DE DE DE DE HE HE HE HE DE IE DE DE DE IE DE HE DE HE DE DE DE IE DE DE DE DE DE HE HE HE DE HE HE HEHEHE HE HE HEHEHE HH HHH HH HHH 
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FILE: Vat FORTRAN A 


‘ SUBROUTINE REQUIRED BY PROBLEM SOLVER, NO FUNCTION IN ANALYSIS 
9300000000000 9 0000980000050 I IIIT II III IIIS III 
INTEGER NEQ 
REAL X,YCNEQ), DYPDYCX) 


RETURN 
END 
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